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Introduction 


This report presents the results of a study of "The Use of Tethers 
for Payload Orbital Transfer" and Is Vflume II of the Semi-Annual Report 
required by the contract. This work was carried out under Modification 5 
of Contract NAS8-33691 originally titled "Investigation of Electrcdynamic 
Stabilization and Control of Long Orbiting Tethers." Dr,, Giuseppe Colombo 
Is Principal Investigator on this contract. The Smithsonian Astrophysical 
Observatory (SAO), studied the dynamic behavior of the tether and the 
Massachusetts Institute of Technology (M.I.T.), under subcontract SVl- 
52006, studied the facilities and systwns required for "The Use of Tethers 
for Payload Orbital Transfer." The results of the M.I.T. study are presented 
in Appendix A of this report. 

The general introduction to the nature and aoplications of the present 
work can be found in the Initial sections of Appendix A. The tody of 
Appendix A contains detailed technical discussions of various tether 
systems. A wiflaerlcal verification of some of the crucial dynamical calcu- 
lations RAde in Appendix A is contained in the SAO work presented in the 
first part of this report. 

Concurrent with this effort, SAO also performed, under Modification 4 
of the same contract, "The Study of Certain Tether Safety Issues" also 
with Or. Col(»nbo as Principal Investigator. The results of that study are 
given in Volunw I of this Semi-Annual Report. 

Tht body of this report has been assembled from the monthly reports 
submlttud under this contract revised and augmented «diere necessary for 
clarity and completeness. This report is intended to stand alone as a 
sim^tiary of the v«rk done on "The Use of Tethers for Payload Orbital 
Transfer." 

The author of this report is Mr. David A. Arnold. The author of 
Appendix A is Dr. Manuel Martinez-Sanchez. 
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2.0 Studj^ Approach to Dynamic Behavior of the Tether 


The SKYHOOK program has bean used to do simulations of two cases 
considered In the M.I.T. study of the use of the tether for payload 
orbital transfer. The purpose of using SKYI^K Is to provide more 
detailed and realistic simulations of the cases considered In the 
theoretical studies done at M.I.T. In particular, there Is the need to 
study oscillation of the system during the various operations. In the 
case of transporting a payload along the wire, the radial velocity Intro- 
duces Coriolis forces that could set up transverse oscillations of the' 
system. These could be a probl«n especially as the payload approaches 
the end of the wire. 

The use of an orbiting tether system by the Shuttle Involves the 
operations of deploys»«^nt and retrieval which may excite oscillation of 
the system. The SKYHOOK program has been used to verify the theoretical 
predictions of the orbits of the Shuttle, tether systnn, and payload 
made by M.I.T. and to determine the extent to which the results are 
Influenced by 11 brat lens of the system. 

3.0 Payload Transport Along the Wire 

One of the cases considered in the M.I.T. study is the transport of a 
payload along the tether from a heavy lower platform to an upper launching 
platform. A simulation has been done using the SKYKiyiK program to study 
the dynamics of the system as the payload moves along the wire. We assume 
that the payload has a means o' grasping the tether and controlling Its 
speed of trans[ Since the net force due to the gradient of the gravi- 
tational and centrlpital forces Is away from the center of gravity of the 
syst«n, the payload will have to be dissipating energy for most of the trip 
from a heavy lower platform to the upper launching platform. 
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The simulation which has been done integrates the motion of three 
mass points - the base platform, the payload, and the upper launching 
platform. The mass of the tether is neglected. A constant transport 
speed of 10 m/sec has been used in this first simulation. The mass of the 
base platform is 300 tons, the payload is 30 tons, and the launching plat- 
form at the top is 10 tons. The altitude of the base is 300 km and the wire 
is 250 km long. The diameter of the wire was set to 2 ram, which is in fact 
not sufficient to withstand the tension load. The only effect of this 
assumption in the simulation is to make the wire more elastic than it would 
be with a thicker tether. For simplicity the integration has been started 
with the payload 1 km from the bottom moving at 10 m/sec. The startup 
phase has been neglected. The radial velocity results in coriolis forces 
that push the payload to the rear. Runs have been done for 100, 500, 2400, 
and 24000 seconds in order to approach the problem gradually in anticipation 
of possible Instabilities. 

The climbing of the payload along the wire has been simulated in the 
SKYHOOK program by making the natural length of each of the two wire segments 
a function of time. The Initial lengths of the lower and upper segments of 
wire are chosen in such a way that they will be 1 and 249 km in length 
respectively when the system is stretched to equilibrium. At later times 
the length of the lower segments is computed as 1^ vt and the length of 
the upper segment is 1] - vt, where v is the velocity of the payload and 
1^ and I 2 are the initial lengths. 

In the simulation, not too much happens in the first 100 seconds. 

Figure 1 shows the results for the first 500 secc^s. Part a) is the in- 
plane motion vs. time, part b) 1s the tension and part c) is the radial vs. 
in-plane configuration at successive time intervals. Mass 1 is the lower 
platform. Mass 2 the upper platform, and Mass 3 the moving payload. Coriolis 
forces »^esult in a displacement of the payload to the rear (positive in-plane 
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Figure 1 Motion during the first 500 seconds with a 30 ton payload 
climbing at 10 ^Vsec starting 1 km from the bottom of a 
250 km wire. Part a) is the in-plane displacenent of each 
mass vs. time, part b) is the tension and part c) (next 
page) is the radial vs. in»plane configuration at 5 second 
intervals. 
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Figure 1(c). Motion during the first 500 seconds with a 30 ton payload 
climbing at 10 ro/sec starting 1 km from the bottom of a 
250 km wire- radial vs. in-plane configuration at 5 second 
intervals. 
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displacement of Mass 3 In Figure la). The upper platform Is generally moving 
forward and closer to the lower platform during this time period. In Figure 
1b we see that the tension In the lower segment Is Initially lower than that 
in the upper segment. This Is because th# center of gravity Is Initially 
about 7.44 km from the lower platform. The tension Is greatest at the 
center of gravity In equilibrium. The payload will Initially have to expend 
energy to get to the center of gravity and will then coast the rest of the 
way up. As the payload moves up, the center of gravity will shift upward 
and the payload will be at the center of gravity 8.06 *'m from the 
lower platform. The tension In the lower wire segment oscillates with de- 
creasing frequency as the payload moves up the wire. The natural period 
for a 30 ton mass at the end of a 1 km Kevlar wire 2 ram In diameter Is 
about 73 seconds. At 5 km the period Is about 164 seconds. These numbers 
agree roughly with the periods seen in the plot. In 500 seconds the payload 
moves from 1 km to '.m from the lower platform. Figure Ic shews the in- 
piane vs. radial configuration at 5 second intervals plotted at equal scale 
in the two axes. We see the slight bending of the wire to the right as a 
result of corlolis forces. 

Figure ?. shows the behavior during the first 2400 seconds plotted at 
25 second intervals. Part a) is the In-plane vs. time, part b) Is the ten- 
sion vs, time, and parts c) and d) are the 1n-niane vs. radial con"' igurat Ions 
at successive time intervals. In part a) we see that the upper bvjSs which 
had been nroving forward for the first 1000 seconds has imved to the rear 
and Is almost in line with the payload climbing the wire. In part b) we 
see that the tension Is now greatest In the lower section since t^»e payload 
passes the center of mass of the system at about 700 seconds. The frequency 
of the tension oscillations Is continuing to decrease as the length of the 
lower section of wire increases. Part c) shows the radial vs. in-plane 
configurations at 25 second Intervals plotted at equal scale In both axes. 
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Figure 2. Motion during the first 2400 seconds of a payload climbing 

the tether. Part a) is the in-plane, part b) is the tension, 
part c) Is the in-plane vs. radial at 25 second Intervals 
plotted at equal scale in both axes, and part d) is the in- 
plane vs. radial with the In-plane axis expanded. 
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Figure 2(d). Motion during the first 2400 seconds of a payload 
clirjing the tether - in-p’'ne vs. radial behavior 
with the in-plane axis expanded. 



Part d) shows the same thing with the In-plane axis ex|:>>;'ided to show the 
in-plane moveinent with better resolution. The features described in part 
a) can be seen in this plot, especially the swinging of the top mass to 
the right near the end of the plot. 

Figure 3 shiMs the behavior for the first 24000 seconds plotted at 
250 second intervals. Part a) is the in-plane, part b) the tension, and 
parts c) and d) show the tn-plane vs. radial configuration. The period 
for in-plane i^endulum oscillations of a tethered system is the orbital 
period divided by the square root of 3. For a 300 km orbit, the orbital 
period is about 5430 seconds and the in-plane period is about 3135 seconds. 
In 24000 seconds we would expect about 7.6 cycles. This seems to agree 
roughly with the results seen in part a). The in-plane period is independ- 
ent of length, so we do not see a change in period with time. In addition 
to the pendulum motion of the system as a whole we also see transverse 
oscillations of the payload on the wire and oscil lotions of the upper 
platform with respect to the payload. When the payload is close to the 
lower platform, the period of transverse oscillations of the payload is 
short, and the period of oscillation of the upper platform is the period 
for oscillations of the system as a v/hole. As the payload climbs, the 
period for transverse oscillations of the payload lengthens, and the 
period of oscillation of the upper platform shortens. Part b) shows the 
tension vs. time. The tension in the lower segment continues to increase 
with time as the payload climbs the wire. The lengthening of the wire 
seen In parts c) and d) is the result of using too sm»ll a wire diameter. 
The cnnputer run halted with the diagnostic that the stepsize was too 
small as the length of the upper segment approached zero. The last output 
point iras at 23750 seconds and the run ended at 23881 seconds. The last 
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Figure 3 Motion during the first 24000 seconds of a payload climbing 
the tether. Part a) is the in-plane motion, part b) is the 
tension, part c) is the radial vs. in-plane configuration 
at equal scale, and part d) is the radial vs. in-plane with 
the in- plane axis expanded. 
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Figure 3(d). Motion during the first 24000 seconds of a payload 
climbing the tether - radial vs. in-plane with the 
in-plane axis expanded. The small oscillation of 
the paylnad and subsatellite as the payload climbs 
the wire is clearly seen here. 



130 seconds have been rerun using the last output state vector as the 
initial state vector as the initial state vector of the next run. The 
results are plotted in Figure 4. The payload has been taken as the 
reference point in the plots so that we see the motion of the upper plat- 
form as viewed from the payload. Part a) is the radial component, part b) 
the in-plane, and part c) the in-plane vs. radial configuration. In part 
a) we see that the upper platform goes below the payload during the last 
couple of seconds. To give better resolution, the last 15 seconds have 
been plotted alone in Figure 5. Part a) is the radial, part b) the in- 
plane, and part c) the in-plane vs. radial. Part c) clearly shows the 
upper platform looping around the payload in the last few seconds. It 
is remarkable that the behavior is stable for so long. The payload would 
of course have to decelerate as it reaches the platform. The rate may 
have to be controlled to eliminate oscillations during the approach to 
the launching platform. 

Since the dynamics of the out-of-plane conqjonent is different from 
that of the in-plane, a run has been done with an initial out-of-plane 
displacement for the wire. The platform at the top was moved 3 km and 
the payload 12 meters placing it in a line between the upper and lower 
platforRts. The wire diameter in this run is 7.6 mm which is sufficient 
to withstand the tension load. The results for the first 23750 seconds 
are shown in Figure 6. The time required to reach the upper platform at 
10 m/sec is 24824.7 seconds. Part a) of the Figure shows the in-plane, 
part b) the out-of-plane, pert c) the in-plane vs. radial, and part d) 
the out-of-plane vs. redial. The in-plane is similar to the results seen 
before without the out-cf-plane displacement. The out-of -plane behavior 
is very regular and does not show the transverse oscillation induced by 
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Last 130 seconds of a payload climbing the tether to an 
upper launching platform. The motion of the upper platform 
is shown relative to che payload climbing the wire. Part a) 
is the radial vs. time, part b) is the in-plane, and part c) 
(next page) is the in-plane vs. radial. 
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Figure 4(c). Last 130 seconds of a payload climbing the tethers to 
an upper launching platfom. The motion of the upper 
platform is shown relative to the payload climbing 
the wire - radial vs. in-plane behavior. 
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(a) (b) 


Figure 5. Last 15 seconds of a payload climbing the tether to an 

upper launching platform. The motion of the upper platform 
Is shoinr) relative to the payload climbing the wire. Part a) 
is the radial behavior vs. time, part b) Is the In-plane 
behavior vs. time and part c) (next page) is the in-plane 
vs. radial behavior. 
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Figure 5(c). Last 15 seconds of a payload climbing the tether to an 

upper launching platform. The motion of the upper platform 
is shown relative to the payload climbing the wire - in- 
plane vs. radial behavior. 
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radial, and part d) the out-of-plane vs. radial. 
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Figure 6(d) 
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CO -i oils forces In the in-pUne component. Transverse oscillations could 
be produced In the out-of-plane component by having the payload out of 
line along the wire. The run was continued another 1070 seconds with out- 
put every 10 seconds. At the last output point which Is 4.7 seconds from 
the end of the ascent, the platform Is 20 meters above the payload In the 
radial direction, and 42 nieters displaced In the in-plane direction. The 
behavior seems essentially the some as In the previous run with no out-of- 
plane displacement. 

In conclusion, the simulations of payload transport along the wire 
using the SKYHOOK program Indicate that the process is quite stable. The 
radial motion along the wlr^ introduces corlolis forces that produce trans- 
verse oscillations in the in-plane, but not the out-of-plane direction. In 
the case of a heavy payload approaching the end of the wire at high velocity, 
unstable behavior would result in the last few seconds. A slowdown phase Is 
obviously required. Additional simulations would be needed to develop an 
appropriate procedure and determine if the rate needs special control at the 
end to prevent the buildup of oscillations duimg the final approach. 

4.0 Simulation of a Payload Launch Using an Orbiting Tether Facility 

The M. I.T. section of the Final Report for this contract (Appendix A) 
describes on page 29 a munerical example nf the Shuttle launching a payload 
using an orbiting tether facility. The Shuttle docks with the tether plat- 
form, transfers the payload, deploys the tether, neleases the payload, re- 
trieves the tether part way such that when the tether system is released 
its center of gravity will be back at its original altitude, and then un- 
docks from the tether system at apogee. After undocking the tether system 
continues the retrieval to the original state. The M. I.T. calculations 
assume the center of gravity of the system remains constant during reeling 
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processes end neglects the effect of llbretlons that are generated during 
the reeling. Simulations have been done using the SKYHOOK program tu see 
the Influence of these effects neglected In the thMretIcal calculations. 

The deployment subroutine In the SKYHOOK program uses equations and 
parameters given on page 9 of HASA TM-X-64963, “A tether tension control 
law for tethered subsatel 1 1 tes deployed along local vertical." The com- 
manded length Is g1v»n by the table on page 10 of the report. The para- 
meters on page 9 are computed for a specific subsatel llte mass and orbital 
angular velocity, and the table of coraraanderi lengths Is for a specific 
tether length and deployment sequence. The table of commanded lengths has 
the undesirable quality of being discontinuous. In order to run the case 
described above the deployment subroutine has been rewritten In a more 
^neral form ^ttemed after the retrieval subroutine. The parameters 
have been rewritten to use the actual masses and orbital angular velocity. 
Instead of using a table, the coinmanded length Is computed as a fraction 
f of the actual length. For retrieval a value of .93 for f gives a slow 
stable retrieval. For deployment, f Is greater than unity. It should be 
possible for f to be substantially greater than unity since deployment is 
an inherently stable operation In contrast to retrieval which must be done 
carefully In order to make sure there is no residual angular monwntum that 
will cause the subsatellite to wrap around the Shuttle during the final 
stages of the approach to the Shuttle. 

The SKYHOOK program has an input paranwter the ejection velocity to 
be used on deployment. This ejection velocity may be large without Intro- 
ducing Instabilities. In this way it Is possible to quickly arrive at a 
sufficient distance from the Shuttle to obtain an adequate gravity gradient 
force for maintaining tension during the rest of the deployment. During 
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the 1n1. n phase the kinetic energy of ejection can be used to maintain 
tension. This tension will eventually use up the initial kinetic energy, 
but by then there is sufficient gravity gradient to continue the deployment 
under positive tension. The tension control law uses the reel motor to 
simulate a visco-elastic tether tuned to the H brat ion frequeficy of the 
tether system. The viscous part of the control law provides the tension 
needed during the initial phases of the deployment. Test runs have been 
done with different ejection velocities to determine a su1tal.«e value for 
running the simulations of the launcn sequence studied at M.I.T. One un- 
desirable aspect of the dynamics is that the control law ends up slightly 
retrieving the system after the Initial kinetic energy has been exhausted. 

In one test run, a wire mass point had just been deployed and the slight 
retrieval caused the wire length to fall slightly below the natural /"ength 
of the wire segment. Since the program is not set up to eliminate mass 
points during deployment, there was loss of tens^on, and the tension control 
law was unable to operate properly. In lieu of pursuing a solution to this 
problem, which would be beyond the scope of this study, the M.I.T. case has 
been run without tether dynamics, integrating only the motion of the two 
end trasses. More study of the deployment process is necessary to useful 
optirriize the process. 

In order to run the deployment, two sets of initial conditions need 
to be computed. The program uses only the state vector for the Shuttle 
initially. The DUMBBELL program is set up to compute initial conditions 
for tv«) or more masses. By making some changes in the program to avoid 
singularities, it was possible to run the program with zero tether length 
to get the initial conditions for the Shuttle. The parameters of the 
system when fully deployed must be given for the other masses. Appendix A 
gives equations for co^uting the parameters of the system at each stage 
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of the operation. Since the SKYHOOK runs have no tether nuiss Included, 
the parameters had to be recomputed with Mj set to zero. The equation 
for 1/L on page 27 of Appendix A is singular for W.|. equal to zero. The 
equation has been rederived without My to get a non-singular expression. 

The first parameter needed to compute Initial conditions for the deployment 
phase Is the value of x, which becomes 21.18 km with no tether mass. This 
places the Shuttle at 378.82 km after the deployront Is completed with the 
upper mass at 478.82 km. A tether diameter of .5 cm is sufficient to with- 
stand the tension load, assuming a break strength of 2.7 x 10*'’ dynes/cm^ 
and a safety factor of 4. Equilibrium parameters for this phase have been 
computed using the DUMBBELL program and used to do the deployment run with 
SKYHOOK. 

The deployment run has been done using an ejection velocity of 5 meters/ 
second. Figure 7 shows the results during the first 2000 seconds at 100 
second Intervals. Part a) 1s the tension vs. time, and part b) Is the radial 
vs. in-plane configuration. Figure 7b uses a new plotting package recently 
developed In which the direction of motion has been reconciled with the 
order In which the configurations are plotted. Successive configurations 
have always been plotted to the right, but in the previous plotting package, 
the Shuttle motion was to the left. For this and all future plots, the 
direction of motion of the Shuttle Is to the right. This change was Imple- 
mented by reversing the sign of the horizontal (in-plane) component of each 
individual configuration. This is equivalent to looking at the orbit from 
the other side so that the direction of motion Is reversed. In part a) the 
tension is Initially high because of the damping term In the control law. 

The control law halts tht outward motion of the subsatellite after about 
500 seconds and there is a slight retrieval during the next few hundred 
seconds as seen In Figure 7b. The deployment then resun^s again. 
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One of the parameters of particular interest in this case being studied 
is the orbital altitude. This information is not contained in the standard 
SKYHOOK output. It can U obtained from the state vector printed at each 
output point. For convenience in restarting runs at a particular output 
point, a special version of subroutine DMPZ is used which writes the state 
vectors on a separate output file. A small program then reads this file to 
find the time requested and formats the state vector for input to a new SKY- 
HOOK run. This formatting program has been modified for this study to also 
compute the radius vector and the magnitude of the velocity from the state 
vector. The altitude is computed by subtracting the earth radius, and then 
plotted along with the velocity using the printer page as a graph. 

Figure 8 shows a condensed plot of the Shuttle altitude HI, payload 
altitude H2, Shuttle velocity VI, and payload velocity V2 during the deploy- 
ment. During the first 1300 seconds which is roughly one quarter of an 
orbit, the Shuttle altitude increases from its initial value of 400.00 km 
to about 400.75 km. The altitude of the Shuttle should, of course, decrease 
during deployment. The initial increase in altitude is the result ;tf a 
slight eccentricity in the orbit introduced by the ejection velocity of 5 
m/sec. This gives the center of mass of the system a radial velocity of 
about 1.06 m/sec which should result in an altitude variation of about .94 
km. Figure 8a shows the orbital eccentricity du»"ing the first revolution 
superimposed on the decrease in altitude resulting from the deployment. 

This eccentricity complicates the interpretation of the results. It could 
be eliminated by giving the Shuttle the reaction velocity that it would 
actually acquire during ejection of a payloaa. 

Figure 9 shows the in-plane vs. radial behavior for the full deployment 
run. The deployment is ccmipleted at about 1800 seconds. In the SKYHOOK 
program, the tension at the Shuttle is computed from the control law during 
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Figure 8. Altitude H (cm) and velocity V (cm/sec) of the Shuttle (mass 
1) and subsatellite (mass 2) plotted at 500 second intervals 
during the deplo 3 mient phase. The period from 18000 to 25000 
seconds is a steady state integration after completion of the 
deployment, a) Shuttle altitude vs. time, b) subsatellite 
altitude vs. time; c) (next page) shuttle velocity vs. time, 
d) subsatellite velocity vs. time. 
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Figure 8. (Cent.) c) Shuttle velocity vs. time, d) subsetellite velocity vs. time 
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deployment. When th: tether reaches its full deployed length, the program 
switches to the steady state mode of integration where the tension is com- 
puted from the tether elasticity and damping. In this run, the damping 
parameter has been set to the value required for critical damping of the 
longitudinal oscillations of the subsatellite at the end of the tether. 

In this way, the momentum of the subsatellite is arrested without recoil 
at the end of the deployment. Although the tether itself has little internal 
damping, the reel motor could simulate a damper if operated 'inder an appro- 
priate control law. At the end of the deployment, the in-plane displacement 
of the tether is about 15 km to the rear, which is an angle of about 8.6 
degrees. After completion of the deployment, the system librates as seen 
at the end of Figure S. The Vibration could be avoiood by introducing a 
control law that terwinates the deployment with a slow-down ohase where the 
deployment rate is controlled so that the wire returns to the ver^'ical posi- 
tion without overshoot. 

The SKYHOOK program terminates the deployment phase at the first output 
point where the tether length exceeds the natural wire length given on out- 
put. The natural length is then recomputed based on the actual length and 
tension at the output point. In this case, the computed natural length used 
after completion of the deploianent was 100.105 km. At the equilibrium ten- 
sion of .6939 X 10^ dynes, the actual length of the tether is 100.6 km. 

Since the system is librating, the tension varies from about .592 to .836 
X 10^ dynes and the length from 100.54 to 100.71 km. The altitude of the 
Shuttle varies from 376.2 to 378.8 km and the altitude of the payload from 
476.3 to 479.0 km. The computed altitudes of the Shuttle and payload fully 
deployed were 378.8 and 478.8 respectively. Comparison of the computed and 
actual altitudes is complicated by the fact that the deployed tether length 
is .6 km too long, the orbit has an eccentricity causing an altitude fluc- 
tuation of about .94 km, and the system is librating with an amplitude that 
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can cause an altitude fluctuation at the ends of about one km. The maximum 
altitude of the system is about equal to the theoretically computed value, 
but the average altitude seems to be on the order of one km lower. It 
might be useful to do a more careful analysis of deployment, retrieval, 
and librations to study possible interactions with the orbital dynamics of 
the center of mass. The output from the SKYHOOK program contains the in- 
formation necessary to compute the work done by the reel motor, the gravi- 
tational potential, the kinetic energy of the center of mass, and about the 
center of mass. The orbital angular momentum can also be studied. 

In order to see the effects of li brat ion, the rest of the study is 
divided into t*<»o case*>. In the first case, the payload was released during 
the forward swing of the tether at the point where the tether is vertical 
and has its maximum forward velocity. In the second case the payload is 
released on a backv«rd swing. The orbit of the payload after release 
requires no marwrlcal Integration and can be calculated from the state 
vector at release. The orbital parameters of interest are the semi -major 
axis a and the eccentricity e. The semi -major axis is given by 

a = l/(2/r - v2/GM) 
and the eccentricity is given by 

e = /I - r^v^/GMa 

where Vj, is the tangential velocity. The program for plotting the radius 
vector r and the magnitude of the velocity v has been modified to compute 
the tangential velocity from the state vector and calculate a and e at the 
time requested on Input. The apogee and perigee are giver, by a + ae and 
a - ae. 

The state vector at 18800 seconds has been used to calculate the 
orbit of the payload released on a forward swing. The payload goes into 
an orbit with a perigee of 478.4 km and a apogee of 1075. C km. For the 



backward swing the state vector at 20400 seconds has been used. The payload 
orbit in this case has a perigee of 476.7 km and an apogee of 896. T km. The 
apogee is almost 180 km higher when the payload is released on the forward 
swi ng . 

The state vectors at 18800 seconds and at 20400 seconds have been used 
as the initial conditions for the second stage of the operation which is 
retrieval of the subsatellite until the center of gravity of the tether 
system is at the original altitude of 400 km at apogee. With the tether 
mass included, the tether should be retrieved to a length of 50.56 km. 
Without the tether mass, using the formula 

l/L = (M;_p * Mup) (Msp * 

the system should be retrieved to a length of 47.068 km (the terms in 
the equation are as defined on page 24 of Appendix A. Figure 10 shows 
the results of two retrieval runs. Parts a) and b} are the tension and 
in-plane vs. radial plots after release on the forward swin?, Parts c) 
and d) are the tension and in-plane vs. radial plots after release on a 
backward swing. The case for the forward swing was run for 7400 seconds 
until the tether was retrieved to a length of 39.17 km. The case for the 
backward swing was run for 7600 seconds to a tether length of 39.38 km. 
Interpolating in the plots of tether length vs. time to obtain the point 
where the tether is 47.068 km long gives 5970 seconds for the forward case 
and 6028 seconds for the backward case. In Figures 10b and lOd we see 
that the initial librations have been damped out and the tether is being 
retrieved at a steady angle which brings the subratellite slightly ahead 
of the Shuttle. An appropriate control law could return the tether to 
the vertical before ending the retrieval if this were desirable. 
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The SKYHOOK runs have been done with output every 100 seconds. In 
order to obtain the state vector where the tether length Is 47.068 km it 
is necessary to interpolate In the output. The program described earlier 
for reading the output state vectors and formatting them for Input has been 
modified to Interpolate between the output points. As a check, the tether 
length Is also computed for the Interpolated state vector. The interpolated 
state vectors at these times have been used as input to the third stage of 
the processing which is steady state integration from the end of the retrieval 
to the next apogee passage where the Shuttle undocks with the tether system. 

In order to determine the orbit of the tether system after undocking from the 
Shuttle the program for reading the state vectors has been nradified to read 
the masses of the upper and lower pallets, and compute the state vector for 

< Vi ^ f C3 ¥ t vji VI V* (%• WV V t I V jr a • Iril^ ▼CVWVI Wf I 

to calculate the post release semi -major axis and eccentricity of the orbit 
of the tether system. 

The first runs done in the steady state phase were unsatisfactory because 
of the linear interpolation used to obtain the initial conditions from the 
output of the retrieval phase. Since both the position and velocity are 
rotating vectors, the linear interpolation results in a systematic shortening 
of the magnitude of the radius vector and velocity, which makes the orbits 
too low. The perigee is reduced by approximately 7 times the error in the 
radius vector which was about 10 km in one of the cases. The interpolation 
has been modified by retaining the same linear interpolation for the direction 
of the state vectors but obtaining the magnitude of the vectors by linear 
interpolation between the magnitudes of the output position and velocity 
vectors. This method should give much better results particularly for a 
circular or low eccentricity orbit. For each of the two cases (payload 
release on the forward and backward swings) the orbit of the center of mass 
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of" the tether system after undocking at apogee has been calculated 
arwilytical ly from the state vector at apogee. The final retrieval run to 
put the tether into its original condition has been omitted since it does 
not app^ear to be essential judging from analysis of the runs up to this 
point. 

An assumption inherent in the theoretical formulas used to calculate 
the state of the system at various stages is that the center of mass of 
the system does not change significantly during deployment and retrieval. 

In a long system, there is a difference between the center of mass and the 
"orbital center" of the system defined as the point where the gravitational 
and centrifugal accelerations are equal. For this case, the orbrtal center 
of the whole system fully deployed is at 399.756 ksn when the center nf mass 
is at 400 km. That is there is a difference of almost 1/4 km between the 
two centers of the system. The angular velocity before deployment in a 
circular orbit at 400 km is .001131402 rad/sec. When deployed with the 
center of mass at 400 km in a circular orbit the angular velocity is 
.001131463 rad/sec. The program for computing the orbit of the center of 
mass has been tested on a short equilibrium run in the fully deployed state. 
The program computed an apogee of 401.46 km and perigee of 400.00 km. The 
distance frcwi the apogee to the orbital center at 399.757 km is 1.70 km 
wh'' :h is 7 times the distance of .243 km between the center of mass and the 
orbital center. Since it is the orbital center of the system that orbits 
the same as a free particle, it has been decided to use the orbital center 
rather than the center of mass as a reference point for studying the motion 
of the center of the system. The position of the orbital center of the 
system r is given by the expression 

r * [i:m^r^./(lm^/rp]^^^ 
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where the r's are measured from the center of the earth. The program has 
been modified to compute the state vector at the orbital center and derive 
the orbital parameters of that state vector. 

In order to study the behavior of the center of the system as a 
function of time, two additional plots giving the altitude and velocity of 
the orbital center have been added to the program that reads the output 
state vectors. Figure 11 shows the plots for the deployment phase of the 
operation. The orbital eccentricity resulting in fluctuations of about .93 
km shows clearly at the beginning of the run There is a decrease in mean 
altitude at the end of the run to about 398.5 km. The decrease in altitude 
is larger than can be attributed to the difference between the center of 
mass and the orbital center. The decrease of 1.5 km is on the order of 
other minor effects and has not been studied to understand the underlying 
reasons. An approach for studying the problem is discussed earlier in this 
report. Plots of the altitude of the center of mass have been done for the 
other phases also. Since the other phases are all less than two orbits it 
is difficult to detemnne a mean altitude from the plots. Orbital elements 
have been computed at the beginning and end of each run for the orbital 
center. In addition, the elements for the tether system and the Shuttle 
have been calculated from the state vector at the time the Shuttle undocks 
at apogee. 

Table 1 gives the apogee, perigee and semi-major axis (SMA) in km for 
the times of interest in all the SKYHOOK runs. Run 1 is the deployment of 
the system out to 100.6 km. Run 2 is the retrieval to 47.068 km after re- 
leasing the payload at 18800 secoi ds on a forward swing of the tether. Run 

3 is the retrieval after release on a backward swing at 20400 seconds. Run 

4 is a steady state integration starting from the state vector after 5970 
seconds of retrieval in run 2. Run 5 is a steady state integration using 


37 . 



TiM (Sec) 


V 

MM A" 

^0.0 

Old*' «• ' , 

400.0 

0 41 « 

1 0 c . 0 

0.1 1 < 

1 400,0 

0 4 ' . • 

ISOC. 0 

d.lD.’.l*- 

;»rc p 

0 . 4 > 0 . 

IU>'0 0 

fl . J “M ' ( ■ 1 

J 11 0 . ft 

0 ) A , ) ; 1 1 ' 1 

«'W0 . 0 

0.1 S > ' • t 

4 . 0 

0.1'*' * ’ 

. : 

M i < . 

4 S t .i . 0 

C 1 VV 

(UCd.O 

w . 4 ' . 

t k) 0 . ' 

k . 4 

V 0 1 . J 

t . 4 : 

-I 

V . ' 

^ V t V ■ . 

V • ■ . 

w ' V ^ .0 

: ' 4 • , 

V 0 0 . P 

(',)»•• 

0 

0. 1 • • ’ . 

itcoc.o 

0 • 

i 1 ^ . K 

1 u • . 0 

) . 1 » . 

1 ' ' : . 

0 . M 

; : . . . c 


: ; i . , 0 


J 1 . V v\ 0 

y 4 

M'.tJO. 0 

0 . • C • . ' 

U '1 d 

C . 1 • 

1 * t V 0 . 0 

0 I < 1 

1"' 'I. t 

C , - ' 

l*Sii 0 

C . V . . 

j • rcc . 0 

0.1 V . . 

u . . 5 

V * ; , 

i * V - . J 

w : • ’ 

1 ' . . . w 

C . • . 

lb. C ' . 0 

k' . ’ 1 ' . , . ■ ’ 

'Oi>.C 

0 . M ' 

1 ' • 1> J 0 

0. • i 

If *100. C 

C 1 » ' • 1 ■ 

jr M 0 

0. ' 

J j , ' . 0 

C . . 

2i . vO.O 

: . t • ; . 

; 11 c j . 0 

C 1 *r. ’ 

23000. <9 

»' 1 : f . . . 

2i;oo.o 

0 . ) « * « ' 

13000.0 

0. JVOt . ' . ■ 

li40C,.0 

0. )»o^' : . 

140CO.O 

0.}> J«. 

S4VC0.0 

0.)OiS«.>i V 

UddO.O 

0.>&74«iU 


Altilurtc (cn) 


OVr.G.N’' 


p,' 


OF FCOt^ quality 




End of Deployment Phase 


Time (Sec) 

w 


Velocity (cm/ sec) 


^ 0.0 

0 . T 6ui'f. 4 1 4' ■ M 

300.0 

0 . ‘ bP ’ • • -1 

icoo.o 

« ' k» <> ' t (1 • 

4 00 . P 

0 , 1 « ' ^ ■ > 

aion.o 

0 1 s r ‘ , f t 

iSOft. 0 

0 . H t > ' ; . t, 

looc.o 

0 16v4< 

3V00.0 

0 , 16' »4i • . . • ^ 

0000.0 

0. * 

4>00.0 

0. >k* .V . • • 1 0 

11000.0 

0.1fre*.'n f 

1000.0 

0. HvU • ( - 

4000.0 

0. »6( It 3 - 

0000.0 

C. »6f ” U- ^ 

0000.0 

C . / ■ ' i t . • ''r 

MOO.O 

O.lvf' e 

0000.0 

c. ?ftc *1 j Ci : M • .'» 

0%«10. 0 

0 . M 1 4 1 ^ ’ r • t* 6 

0000.0 

C.7v4l’-ft 'ui 

0000.0 

0. ■ 6' • ' ti ti. • t 

ioooo. 0 

C . Tie k -C*. . ' ’ • 1 

40500. 0 

0. J : ^ ‘>1 4. M)4 

ilOOO.O 

0, }ttu Mil. 

nooo.o 

0. ’46r M ’*.'11 MJ4 

llOOO. 0 

0 . ?ra 5 - : -o< 

13500.0 

0 . ’ t 2 4M*. J . t 

UOOft.C 

0 . >*>p 'fv- )* , jr •• 6 

41500.0 

0 . pf. 

14000.0 

o.T»t„ k^. n, •, ..-a 

14000.0 

0 . T f f n ; . * 3 ” i. . . % 

10000.0 

6 . T6f fl ( . 4 ; r,. . - , 

10500. 0 

0 . ;t<.«dSuL s ,tm 4 

10000.0 

Q. .n* 

10000.0 

0 . 7 4ft?) U )t M..Ct 

11000.6 

0.7 4444 74 1 ; ' r • kit 

nooo.u 

Q. 74«o; M 4 *06 

10006.0 

0.I444M 4 ^t)U-n»i 

10000.0 

0.7««b > SVSi *t >r 6 

10000.0 

0 . 7tC4 Jk- ’ * 2 '.. M 6 

10006. 6 

0 . 7 4 7 c 3 A 1 * •. 1 • d t 

aoooo.o 

0.7470 Jb> ll.*04 

Oosoo.o 

0.74M0r*M .'IL M14 

11000.0 

0.7470514 « ) :i 06 

11000.0 

0.74447 ) 1 'd U .06 

laooo.e 

0.74<.»47 7fr>t! 04 

aaooo.o 

0.7044)42«PU«0« 

aaooo.o 

Q.I441itOS77L*C6 

nooo.« 

o.m4im)4so4 

aoooo.o 

0.mift«02U:«04 

00100.0 

o.i04»}o;7a«K«o4 

00000.0 

0.l0m01l4dO«04 




Figure 11. Altitude (cm) and velocity (cm/sec) of the orbital center 
during the deployment phase and steady state integration 
just after completion of deployment. 
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Lirw Run 

Time 

Masses 

Apogee 

Perigee 

SMA 

1 

1 

0 

86+23.109 

400.9365 

359.0645 

400.0005 

2 

1 

18800 

86+23.109 

399.5903 

356.9313 

j 18. 2603 

3 

1 

10800 

19.109 

1075.0220 

47B.4437 

776.7329 

4 

1 

20400 

86+23.109 

400.2854 

397.2581 

398.7717 

5 

1 

20400 

19.109 

896.1484 

476,7392 

6G6.4438 

6 

2 

0 

86 + 4 

302.1875 

265.7070 

323.9473 

7 

2 

5970 

06 + 4 

302.5229 

266.0026 

324.3028 

0 

3 

0 

86 + 4 

380.6198 

301.6923 

341.1561 

9 

3 

€206 

86 + 4 

300.9478 

301.8230 

341.3054 

10 

4 

0 

86 + 4 

382.5229 

266.0826 

324.3028 

11 

4 

5100 

86 + 4 

382.5223 

266.0773 

324.2998 

12 

4 

5100 

6 + 4 

404.6969 

379.3752 

3S2.0361 

13 

4 

5100 

80 

380.3441 

251.5480 

315.9460 

14 

5 

0 

86 + 4 

300.9478 

301.8230 

341.3854 

15 

5 

4050 

66+4 

380.8854 

301.8820 

341. 3837 

16 

5 

4050 

6 + 4 

419.8848 

396.5458 

408.2153 

17 

5 

4850 

80 

370.7827 

287. 4952 

333.1355 


Apogee, perigee, and semi-major axis at various stages 
of the launch sequence. The values are for the orbital 
i:enters of the masses listed in the fourth column, namely 
he Shuttle (80 tons), lower pallet (6 tons), upper pallet 
(4 tons), upper pallet plus payload (23.109 tons). Shuttle 
plus lower pallet (86 tons), and payload (19.109 tons). 

Run 1 is the deployment, runs 2 and 3 are partial retrieval 
after release on the forward and backward swings respec- 
tively, and runs 4 and 5 are steady state runs from the end 
of the retrieval in runs 2 and 3 respectively to the next 
apogee passage. 


as input the state vector after 6. '06 seconds of retrieval in run 3. In 
run 1, the lower mass is 86 tons and the upper is 23.109. In the other 
runs, the lower mass is 86 tons and the upper is 4 tons, for runs 4 and 

5 the state vector for the orbital center of the tether system alone has 
been computed at the time of undocking. The mass of the lower pallet is 

6 tons and the upper is 4 tons. Finally, the state vector for the Shuttle 
(80 tons) at the apogee of runs 4 and 5 has been used to get the orbit of 
the Shuttle after undocking from the tether system. The masses listed in 
the Table indicate which configuration is being computed. 

Line 1 in Table 1 is the orbit at the beginning of the deployment. 

The computed eccentricity agrees with the plots of altitude in Figure 11. 
Lines 2 to 5 are the orbital parameters at the time of payload release on 
the forward (18800 sec) and backward (20400 sec) swings. The average semi- 
major axis of lines 2 and 4 is 398.5 km indicating a drop of 1.5 km during 
deployment. This agree* with the results seen in the plots of the orbital 
altitude vs. time. The semi-nwjor axis of the orbital center is about .5 
km higher on the backward swing (line 4) than on the forward swing (line 
2). Lines 6 and 7 give the orbital parameters at the beginning and end 
of the partial retrieval after release on the forward swing at 18800 sec. 
The semi-major axis has increased by .36 km during the retrieval. Lines 
8 and 9 are the corresponding results for retrieval after release on a 
backward swing. The increase in altitude here is .23 km. The semi-major 
axis is about 17.1 km lower in run 2 than in run 3 as a result of the 
'greater er*e*"gy given to the payload by releasing on the forward swing. 

Lines 10 and 11 give the orbital parameters at the beginning and end of 
the steady state run from the end of retrieval to tether system release 
at apogee for the case of payload release on the forward swing. The semi- 
major axis is nearly constant. Line 12 is the subsequent orbit of the 
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tether system after undocking, and line IJ is the final nrhit of the 
Shuttle. Lines 14 through 17 give the corresponding information for the 
case of payload release on a backward swing of the tether. Lines 14 and 
15 show no change in semi-major axis. The average semi-major axis of the 
tether system after undocking obtained from lines 12 and 16 is 400.1 km. 

This is within .1 km of the theoretically calculated altitude of 400 km. 

The orbit of the tether system is 8 km higher than predicted for the case 
of payload release on a forward swing and 8 km lower on the backward swing. 
The tether system orbits are eccentric by 12.7 ano 11.7 km for the forward 
and backward cases respectively. The average perigee of the fihuttle after 
undocking from the tether obtained from lines 13 and 17 is 269.5 km. This 
agrees within .1 km with the theoretical ly calculated value of 269.4 km. 

The final Shuttle perigee is 18 km lower or higher depending on whether the 
payload is released on the forward or backward swing. 

In conclusion, the simulations done with the SKYHOOK program give good 
agreement with the theoretical ly calculated results from the M.I.T. study 
and indicate the order of magnitude of the perturbing effect of librations 
not considered in the theoretical study. The results indicate some altitude 
changes during reeling operations, and fluctuations within the libration 
cycle. These effects are about an order of magnitude smaller than the 
impact of releasing the payload on the forward or backward swing. 
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Introduction 


Rocket propulsion is a well established method for moving payloads in 
space, and has thus far enjoyed a virtual onnopoly in such applications. 

This pre«ninence is likely to continue for the foreaeablc future, but 
with the advent of new and more flexible tools for access to space, such 
as the Space Shuttle, and with the constant struggle for n»re payload per 
unit cost, we ar» likely to witness the development of alternatives and 
supplements that will achieve the aame mission with less mass (and cost) 
expenditure. Planetary gravity assist can be regarded as one such alter- 
native, in practical use today; other concepts, such as geomagnetic pro- 
pulsion, interplanetary ramjets, etc., have been proposed as future develop- 
ments. 

In this report, we explore a relatively sioiple concept for enhancing 
interorbital transfer capabilities. It is well known that in an extended 
orbiting body only certain points (those on the Earth-centered circle 
through the body's orbital center) are in centrifugal-gravitational equi- 
librium. Other points in the body undergo a net resultant force (the 
gravity gradient force), which, for elongated bodies, tends to align them 
along the local vertical. Thus, if a satellite is joined to a larger 
spacecraft in circular orbit by means of a long, lightweight cable (tether), 
its equilibrium position would be directly above or directly below that 
spacecraft, along the local vertical. A certain point (close to the system 
center of mass), intermediate between satellite and spacecraft, would be in 
true orbital equilibrium, while the two and messes would be pulling on the 
tether. If the masses of the two bodies ere m; and mj , the cable length 



2 


Is L and the orbital angular velocity Is Q, the cable tension Is 

D] ni2 

F ■ ; 012 " 7 (reduced mass) 

ID ] T" 02 

which can he recognised as the local weight of a body of mass B12 

(close to the smaller of the tvro masses), times a factor 

This force arises because the linear velocity of the (outer) mass 

(inner) 

Is (greater) than that required for orbital equilibrium at Its 
(smaller) 

location. Therefore, If the cable is cut, the upper body will entur 
an elliptical orbit whose perigee Is at the initial altitude of that 
body In the compound structure; simultaneously, the lower body enters 
another elliptical orbit, this one with Its apogee at the Initial 
altitude of this body. 

We can now imagine the lo%rar body to be either an orbiting Space 
Shuttle or an orbiting space platform, and the upper body to he a 
relatively light satellite; if the tether is long, the satellite will go 
into a high elliptical orbit, while the platform will enter an orbit 
with a slightly lowered perigee. If we want the system to be reuseable, 
some thrust must be applied to the platfom to raise It to its original 
orbit. We will see in what follows that the total impulse required to 
do this is at least equal to that which would be required to place the 
satellite In its high ultimate crblt using thrust instead of the tether. 
However, and this is the crucial point of the system, the platform can 
be raised using high specific Isipulae, low thrust, electric rockets in 
the period between missions, while In each mission the satellite is 
inserted In its orbit in a rapid, qua^ii -impulsive manner by the tether. 
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Alternatively, In the case where the Shuttle Is the lower spacecraft, we 
may choose to reenter after releasing the payload, with no need to reestab- 
lish the orbit. Some operational complications arise In this case, as will 
be discussed in the text, but the maueuver is feasible within certain limits. 
The net effect in either case can be a large overall propellant savings in 
the upper stage, or, equivalently, the ability to transfer significantly 
larger payloads with a given amount of propellant . 

In the present study we have identified and asseesad a variety of 
ways to implement these concepts. Given the limited scope of the effort, 
the study has been restricted to conceptualizing each system, performing 
first order orbital calculations to determine payload gains, and, at least 
for a few of the most attractive systems, carrying out conceptual designs 
that allow estimates of mass and power to be made. In addition, detailed 
tether dynamical studies were made for one case where apparently no prior 
work existed (see Appendix 4 ) , and very limited cost estimations 
were made for some of the systems in order to gain insight on leturn-on- 
Investment times. 

For the purposes of this final report vtt have chosen to present the 
material in what seems to be the most logical order, namely, from the 
simpler and most clearly feasible systems to those whose inq>leiaentation 
offers difficulties, but which, by the same token, offer the greatest 
potential promise in terms of performance. This is not Che order in which 
the work evolved chronologically, and some unevennees of detail may be 
unavoidable as a consequence. 

For Che sake of clarity, we present here a brief description of each 
of the systems studied, with some coments on their salient features. A 
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detailed diacussion of each of thete le to be found in the body of the report : 
Tethers as Shuttle fncilitiee . The aiopleat iaapleaentation of 
the tether concept is when the tether system is pensanently attached to 
the Shuttle and is flown into and out of space in each mission. Clearly, 
this displaces some Shuttle payload » and its utility must be restricted. 

We found, however, that for high energy Shuttle missions, such as polar 
orbits, where payload Is not limited by structural considerations, the 
Shuttle flight envelope can be appreciably extended by a short, rewindable 
on-board tether. For 26 . 5 " orbits, no advantage was found using tethers. 

(b) Space-baaed, low mass tether systems . Another 
promising system is one where the tether and its associated hardware are left 
in space after each reuse. For tether lengths not cxcseding some 100-130 Km, 
depending on payload, the lower mass can be provided by the Shuttle itself, 
docked to the rewinding end of the tether. Payloads are brought up by the 
Shuttle, each attached to its Orbital Transfer Vehicle; they are mounted 
on the tether end, the tether is deployed and the payload is released, after 
which the OTV fires to place the payload in its transfer orbit. The Shuttle 
now enters a lower elliptical orbit, but not low enough to force reentry; 
the tether is partially rewound and released at a condition such that, after 
autonomous completion of rewinding, it is back in its original orbit. The 
Shuttle now reenters. 

For the lengths indicated for this and the previous syst«a, the tether 
mass is fairly small, and winding -unwinding operations, using rate controls 
that have been studied elsewhere (Ref. 11 ), should present little problem. 
Payload Increases of some 13% are indicat id for a Centaur OTV used from LEO 
to GEO with a 100 Km tether. Deep space iselons can alao be significantly 


enhanced . 
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(c) Pl«tform-ba8«d Infraedlate fther syteas . Highsr pcrfo'.iaance 
can be obtained with longer tethers. At this point, however, a lower platform 
more massive than one empty ShutLle becoiaes necessary to prevent reentry after 
release of the payload. Orbiting space stations of the kind envisioned for 
the 1990 's are natural candidates. Now, the platform orbit must be restorcu 
by application of low level, long duration, high specific Impulse thrust. 

This, in turn, establishes fairly high requirements for electrical power on 
the space station, which stay become the factor limiting ' achievable tether 
length. In addition, the tether itself becomes too massive to be conveniently 
rewound sfter each miasion; an alternative concept that was evolved consists 
of a "ferry" or elevator which travels the length of the tether (up to some 
23C~300 Km) to deliver the payload~OTV combination and return. The dynamics 
of this travelling ferry was studied in some detail, and no real problems 
were encountered, although, as in other systems, climbing rate must be care- 
fully controlled. For Centaur transfers from LEO to GEO, a 250 Km tether 

of this sort allows some 38% payload Increase, but requires about 400 Kw of 
electrical power in the space station (for orbit racovery In 14 days) . 

(d) Large-scale tether systeias for LEO-GEO transfer . We also studied 
more arobiclous systems Involving two permanent tethers, one In LEO (radially 
out) and one in GEO (radially in) . The payload-OTV is released by the lower 
tether and a first Impulse is applied by the OTV to enter a Hohmann orbit. 

At Its apogee, a second impulse matches speed with the lover end of the GEO 
tether, and, after capture by It, the payload travels along the tether to 
GEO orbit. By proper choice of parameters, the rendesvous can be attinapted 
again after an integer number of orbits. This system can In principle be 
pushed to the limit where no propulsion Is needed on the peyload , If the 
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initial LEO orbit is equatorial; however, this requires tethere with lengths 
of the order of 1200 Kin in LEO and 10,000 Km in GEO, and of great nuiss. 
Intermediate solutions are possible using nonzero impulses in LEO and GEO; 
for instance, a 430 Km LEO tether (weighing 7.5 times the OTV mass) and a 
5900 Km GEO tether (of mass 10 times that of the OTV) can be combined to 
obtain a factor of 2.8 in payload capacity for a Centaur vehicle. 

Ref. 1.1 Charles C. Rupp, NASA TT1 x-64963 "A tether tension control 
law for tethered subsatellites deploynt ■'long the local 
vertical . " 
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Tethers as Shuttle Facilities 


Use of a tether system as a permanent facility of the Shuttle does 
not appear Justified for missions that fall within the operational envelope 
of the orbiter with its integral (X1S tanka. This ia because, even though 
the tether allows deployment of the payload from a lower Shuttle orbit 
(typically an elliptic one) , the payload cannot be increased due to other 
constraints, such as payload bay structural integrity and c.g. location. 

Tlje only savinps are then in the use of less (MS fuel, but these cannot 
balance the loss of revenue from the payload displaced by the tether itself. 
An example i.s shown in Tabic 2.1: a 47 Km .«ther allows payload to >e placed 

in a 500/500 Km orbit from a Shuttle in a 185/453 Km orbit, with an OMS 
fuel savings of $33,000. However, the mass and length of the tether 
facility displaces payload worth $2.80 M. Similar results are shown for 
a polar orbit. 

There are some possible scenarios where a Shuttle baaed tether could 
be cost-effective. These refer to low Earth orbits high enough (particularly 
for polar orbits) that payload is limited by (MS fuel capacity, including 
extension kits. A trade-off study is next presented to determine how far 
the operating envelope can be extended by a permanent Shuttle tether. 

If the OMS fuel extension kits are not available to the Space Trans- 
portation *=ystem, then the relevant comparison is between the basic Shuttle 
with only the integral OMS tanks and the Shuttle with the on-board tether 
system. The advantages of the tether system are then apparent, resulting 
in a flight envelope comparable to that afforded by the fuel kits. 
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TABLE 2.1 

COST TO LOW ENERf.Y MISSION** 



Space Telescope 

Polar Orbit 


Orbit 500Km/28.8“ 

1000 Kil/97‘’ 

Weight of Pay loan vkgy 

11,000 

3,000 

Length of Payl,jad (m) 

13.1 

9.0 

Diameter of Payload (m) 

A. 26 


Cost to current Shuttle ($M) 

20.20 

23.07 

Cost to Shuttle + Orbiter based 



tether system ($K) 

23.00 

29.8 

Lost revenue from displaced payload ($M) -2.80 

-6.73 

OMS fuel savings ($M) 

(0.033) 

(0,083) 

Benefit of using tether system 

($M) -2.77 

-6.647 




1) Cost per Shuttle flight "= 

$27.3 at ETR 



$46.9 at WTR 



2) Elliptic Shuttle orbit + tether transfer 
perigee altitude = 185 km 
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Calculation of the Shuttle Mission Envelope with an On-Board Tether System . 

We adopt as a prototype mission one where payload is to be delivered 
to a circular orbit by releasing it from a tether which la attached l:o a 
reeling device carried on board the Shuttle, and which is subsequently re- 
trieved by it. The Shuttle itself flies an elliptic orbit with Its apogee 
lower than the payload orbit (by an amount equal to the tether length) . It 
is injected into this orbit by a velocity increment applied by the OMS at 
some point in a standard 185 Km circular parking orbit. After release of 
the payload, the Shuttle enters a new elliptic orbit with perigee lower than 
the original 185 Km, and then de-orbits by application of an additional AV 
(such as to give a theoretical perigee of 0 Km). For calculations, one 
further OMS firing of 30.5 m/sec (106.7 m/sec for WTR launch) is assimed 
for insertion of the loaded Shuttle into the parking orbit, and a 12.8 m/sec 
AV reserve is assumed. 

An example launch sequence (for an intended payload orbit of 600 Kn 
altitude and 28,5” inclination) is shown in Table 2,2. 

Starting from the parking orbit at radius (equal to the perigee 

length of the transfer orbit), the Shuttle enters the transfer orbit of 
apogee by an OMS firing having a AV of 


AV 


insertion 


y_ 

R 


\V 


V 


- 1 


( 1 ) 


The angular velocity at apogee can be expressed as must 

be equal to that in a circular orbit at the final payload radius 
^T’ tether length. This gives 




2R 

E__ 

R.-^ R_ 



( 2 ) 
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which can be rewritten as 


2R^ 



( 3 ) 


This equation is to be solved for the tether length L,j, in terms of 

the parking orbit radius and the final orbital radius. AV. then 

ins. 

follows from Lq . (1), using • 



11 


Table 2.2. Typical sequence of mission events 

STS/on-board tether system 

From tlR to 28.5® 600 km circular 
orbit 


Event delta-V(OMS) 

m/sec 

Shuttle insertion 30.5 

Injection burn 100.4 

Payload release orbit 
Shuttle after P.L. release 
Shuttle daorblt 55.1 

(note: delta~V reserve * 12.8 m/sec) 


Resultant 

h /h , km 
a p 

185/185 

534.3/185 

600/600 

516.2/68.8 

516.2/0 
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The tension T in the tether is constant If its own mans can be Ignored 
relative to le satellite mass. If we further assume that the satellite 
mass is small compared to that of the Shuttle (a conservative assumption), 
and that the tension at perigee can be estimated as if the orbit were 

circular with R «• R , we can write (Ref. 2.1) 

P 


\ R 


V 


X 


1"+X \ 1+A 


v'h 


ere is the payload mass and \ ■ L.J./R . 


(4) 


For a Kevlar Aramid cable of 0.9 mm diameter the mass is m ■ 0.59 Kg/Km, 

* 

and the minimum break strength is T_ • 90 Kg « 882 Nt . Allowing a 

or • 

safetv factor f (f *=3 was used in calculations), the number of strands 

s s 

T.f 

g 

is , and the tether mass is 

^Br 


T.f m 

M_ - L. (5) 

Br 

The mass of the reeling v..id other devices is expected to be proportional 
to the tether mass. Based on Ref. (2.1), we estimate for the total tether 
system 


\s - ^-7 M,r 


( 6 ) 


Notice that, to this approximation, is proportional to M^, the 


Ref. 2.1 NASA CR-132780, Appendix D, p. 222. 


Note that somewhat different values of the working stress have been used 
in other sections of this report. The value used In each case has been 
clearly identified, however. For a discussion of the uncertainties in an 
estimate of this parameter, see Appendix 1. 
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payload mass; the quantity 


"Vs “ s: 


Ts 


( 7 ) 


is therefore a function of the orbital parameters, but not of 


The mass of the orbiter at ME burnout is 


M 

o 


M 


oe 


M, 




+ M 


( 8 ) 


where M is the empty mass and M the mass of CWS fuel carried (either 
oe p 

that in tlie Integral OMS, or including an integer number of (MS kits as 
well) . The fuel required for Insertion into the parking orbit plus 
injection into the transfer orbit (total velocity increment ■ h\’i) 
is M^(l-ui) where 


Ml 


^-AV,/gI 


sp 


( 9 ) 


The mass after releasing the payload is then M^M i ~ , and after the 

deorbiting burn (AVj), the mass is 


(M^Ui - M^) P: , with 


Mz 


g-AVz/gl 


sp 


( 10 ) 


This mass is to be equated to » since, by assumption, all 

the OMS fuel has been used up: 


\ + "ts * "p>“‘ - "l'"' ■ "oe* "ts 
T s “ "l “ts * 


and using M, 



This equation allows calculation of the maximum payload capability 
of a Shuttle-tether system combination for a given payload orbit and a 
given OMS option. An additional limitation was used, namely, that ♦‘he 
sum of the pavload , tether system and CMS kits should not exceed 29500 Kg, 
the full load capability of the Orbiter. 

The results of these calculations are displayed in Fig. 2.1 (for ETR 
launch into a 28.5° inclination orbit) and Fig. 2.2 (for WTR launch into 
a 104° orbit). The basic Shuttle envelopes shown for reference were 
calculated from the same basic equations, modified to allow variable 
parking orbit” altitudes and no tether. 

Fig. 2.1 shows that the combination Integral -OMS-plus-tether has 
substant iall V more payload capacity than the basic Integral OMS Shuttle 
for a 28.5' orbit; it allows, for instance, 21 Ton payloads to a 700 Km 
orbit without any extension kit. The same is true when OMS kits are added 
to both systems (without and with tether) . When we consider the envelope 
of the curves for the tether system with varying numbers of extension fuel 
kits, we find Chat it does not exceed the corresponding envelope without 
tethers. Thus, if extension kits were available, the usefulness of the 
on-board tether vjould be marginal . 

For the polar orbit case. Fig. 2.2 shows the same, and even more 
pronounced, gains in payload-altitude performance for polar orbits. In 
fact, even the envelope is now extended; l.e., certain missions which 
are simply not accessible to the Shuttle with any number of extension 
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kits now become feasible. This is the case for orbits betwaec 500 and 700 
Km, for which payload extension of some 2 tons become possible using the 
integral QMS tanks plus a tether system. An example of the use of this 
extended capacity would be the possibility of placing the 11 ton Space 
Telescope in a 600 Km polar orbit. Similarly, i^lle no payload can be 
delivered by the Shuttle to orbits higher than 960 Km (even with two extra 
fuel kits), the tether system with one single kit allows some 4.4 tons to 
be placed in 1000 Km orbit (2 tons with the tether and no kit) 




Circular orMfc altitude (m) 










Circular orbit altitude (rn) 



18 


3 • Space-Based, Low Mass Tether System . 

Prelim, tnary CoosiderationB . 

The concept of leaving the tether in orbit for reuse was introduced 
early in our study, and validated by simple orbital dynamics calculations, 
which showed marked increases in payload capability both for orbit trans- 
fers and for deep space missions. In these initial calculations, the 
reaction mass attached to the tether base was assumed large for simplicity, 
and no account was taken ot tether mass (although the tether length was 
restricted to less than 400 Km to keep its mass within reasonable limits, 
see Appendix 1) . 

Table 3.1 shows some results of these calculations, assuming a two- 
stage lUS vehicle is attached to the payload and used for the Initial and 
circularization firings in a LEO-GEO transfer. The tether is attached 
to a massive LEO base. Payload increases of roughly 20% per 100 Km of 
tether are predicted. 

Similar results for deep space missions are shown in Table 3.2, this 
time in conjunction with a Centaur vehicle. The value of c^ “ 80 Km^/sec^ 
is typical of direct Galileo orbits, and, as shown, an 8% payload increase 
is predicted per 100 Km tether length. For other excess hyperbolic 
velocities, the results are given in Fig. 3.1. 

3 . i System description . 

While these calculations clearly show the desirability of such tether 
systems, they ignore the complications due to the finite mass of the lower 
platform. In particular, these can be importani if this platform is simply 
the Space Shuttle, plus possibly a lighter station at the lower tether end. 
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Tabl* 3.1 

rAYi.oAD Btm riT r oj< GtocwnmoKcxjs orbit transfer* 


't)r ' 
(laJ 


0 


ion 


200 

300 


400 


PayloaJ Weight 
(kg) 

24Gb 

3122 

3075 

4320 

5100 


Payload incrcbi' 
(V) 

le 

39 

C3 

93 


C&lctildtion coiidititjns : 

1. SnU'i'TU: + Two stage lOG 


2 . 


3. 


Stage 
Isr(Bec) 
f Btru. 

VCT pro}), (kg) 
ParJeing c^rbit: 
'lather cystem dock 


1 2 
291.9 289.7 

.946 .933 

9707 (2722) 

300/300 kin 

with ehuttlc in parking orbit. 
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Tabu 3.2 

p;»Yj.:Ari r.::wEr’iT roi< sou^iR bxmoratjon** 


C3 

Tcliier length 
(km) 

Injected :iiasr 
(ky) 

Inert'd sc 
(t) 

IS 

0 

76‘>3 



inn 

8253 

7.2 


200 

8057 

15 


.300 

9.511 

23.6 


400 

10219 

32.8 

CO 

0 

224«^ 



100 

2413 

7.4 


200 

2589 

15.2 


303 

2771 

25.3 


400 

2963 

31.9 


Calculrttiun conditions: 

1. SHUTTLE + CE!1TAUP. 

Isp * 444 sec 

WT of prope llant «= } 3608 kg 
Dry vrr •- 1827 kg 

2. Parking orbit: 300/30*i )» 

3. Tother system dock with shuttle in parking orbit. 

4. C3-V^-fJ:S 

r 




enhE.nce9 3.atinch 
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This combination is quite attractive in that it reduces to a mlnlmuin the 
needs for elaborate and costly space platforms, and is therefore amenable 
to an early implementation. In the remainder of this section we describe 
and analyze a reuseabie, low mass tether system for use in conjunction 
witli the present Space Transportation System. 

3 . 2 Space-based, low mass tether systems for orbital transfer assist . 

The core of this system is a pair of relatively light space platforms 
connected by a tether of up to about 100/200 Km length. The lower platform 
can be quite similar to the pallets used as Airbonve Support Equipment 
(ASE) for mating the lUS rocket vehicle to the Shuttle payload bay. It 
would be designed to house the wind/unwind mechanism and control r, to house 
the fully wound tether during initial launch and between missiors, end to 
dock and interfac with the Shuttle for subsequent missions. 

Th3 second, or upper platform, has as its mission to receive the OTV/ 
payload package from each loaded Shuttle after its docking with the lower 
platform, to hold this package during tether unwinding, and to release it 
after stabilization at the fully extended position. Due to the low gravity 
gradient forces involved in this system, this upper pallet can be consider- 
ably lighter than the aforementioned ASE. 
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After releasing the OTV, the tether would be rewound in stages as 
discussed below, and the whole system would be left in orbit for reuse. 
Since its total mass is of the order of 15 Tons, the systiMB can be de- 
livered by one single initial Shuttle flight. Its reuseability is in 
principle only limited by tether wear. 

Following is a step-by-step description of the typical mission 
for this system; 


Stage 0 ; 


Stage 1; 


Stage 2 ; 


Stage 3 ; 


Shuttle flight delivers tether system to orbit (between 
300-400 Km) . System consists of a lower pallet, designed 
to dock with subsequent Shuttles and to wlnd-unwind the 
tether, a length of tether (100-200 Km, depending on pay- 
load), and an upper pallet, or teleoperator, designed to 
hold the OTV and payload. 

Later, another Shuttle flight docks with tether system. 

OTV + payload is transferred to upper pallet. Tether is 
unwound slowly, at controlled rate. After stabilization, 
OTV Is released. 

OTV fires, places payload on transfer ellipse. At GEO, 

OTV circularizes. 

Shuttle, docked to pallet and with extended tether enters 
an elliptic orbit, with perigee above sensible atmosphere. 
While in this configuration, tether is partially rewound, 
until its c.g. coincides at apogee with original c.g. 


altitude. 
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Stage A ; Shuttle releases (at one apogee passage) tether system, which 

I 

stays in original circular . t. After release, pallet 
completes tether rewinding. Shuwtle itself goes Into slightly 
modified elliptic orbit, from whlcK it reenters as desired. 
Tether system is ready for reuse. 


^ ^ Performance analysis . 

Nomenclature : 


L 

M 

up 


^SH 
\eo “ 


G 

G' 



M • sum of these 
tot 


*TS 


tether length 
mass of OTV 4- payload 
mass of upper pallet 
mass of tether 
mass of lower pallet 
mass of lower platform (Shuttle)' 
radius of orbit for autonomous tether system. After Shuttle 
docking and tether deplo 3 rment, orbit radius for 

the overall c.g. 

Overall c.g. before OTV release 
Overall c.g. after OTV release 
distance from Shuttle to G 

distance from Shuttle to G' before partial rewinding 

distance from Shuttle to G' after partial rewinding 

distance from Shuttle to tether system c.g. after partial rewinding 
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Before payload release (but after tether deployment), the overal c.g. 
(C) is at a distance x from the Shuttle-mated lower pallet: 



L 


/»L 



«TOT- "SH-^ “l *>^*%*\ 

P 

After separates, the new c.g. (C) is at 

Point G * now enters en elliptic orbit 
with apogee 

V \eo‘ 

and apogee velocity 


V 

a 





Using 


,, 2R 

,2 « E — 


a R R + R 
a a p 


, and expanding 


8t 


to 1 order in (x-x’)/Rj^^q, we find from the above 


“p ' *LEO - ’ 
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Next, we slowly rewind, while G* stays in the aaoe elliptic orbit. 
When the new tether length is £., the distance between G* and the Shuttle 
is 


up T 



T 


where Rj, ” i./L 

\p‘ ^ ^ 


U ' ;Lt5 so that Rj, •« \p - 


jm^A 


Ai.ao, the distance to the c.g. of the 

tether system alone is 


M + Kj,/2 


. -HE 


*TS ■ 


The orbital eccentricity € *■ 


R - R , 

^2 E ~ e X - X* 


R + R 
a p 


forces in-plane 


oscillations of the tether at the orbital frequency and with amplitude e. 

It can be shoi?n readily that their cffect\<^f second order in — ^ , 

— - R^gQ 

and will not be included in this analysis (although they should be assessed 
in a more careful study) . 

The forward speed of G^^ at apogee is therefore 


V 4- (x__- x') iii 
GTS a tS K 
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We want to stop rewinding and release the OTV at apogee when 
coincides with the orbital speed at the locatla i of i at 

*CTS' ^S- *' ■ \e0* V- S' - * + 

This leads to the condition 

-\/ w Ss- « * VS” 

or, after expansion and simplification. 

In words, the c.g. of the tethei system alone must be made to coincide 
with the original overall c.g. If thl^t condition is satisfied at the Instant 
the Shuttle detaches from the partially rewound tether system, the latter 
(its c.g.) will remain in the original circular orbit. Final rewinding after 
this time will not affect this result, and so the fully retracted tether 
system Is ready for reuse. 

Using the formulas derived for x, x*, x,j,g and x* , we can now calculate 
the required partial rewinding length: 

M_ ft > 0 

.M +r^~ M + ~r M2 + iL./2 

* up 2 L _ up 2 L . ^ ^ up T _ _u£ 

\ Mt ^OT" \ 

or, after simplification 

E -up , ..^up/, "tS 

L «I V\' «tOT “t “sh 

”ts ■ "t M ^ (tether system mass) . 


where 
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After relcar-ing the tether system, the Shuttle itself enters a new 
elliptic orbit with apogee at 


R 


- R^, - x' - Rj^gQ 


- X + x' “ X 


with apogee velocity 



Hence its new perigee is at 


"p.SH ■ \e0 - 

This may in some cases be actually higher than the altitude of the Shuttle 
It) the first perigee passage after payload release, but before any rewinding. 
Ituis, one should also check this altitude. 


h R 

M1N,SH. p,G 


“ 4.E0 ■ ^ 


. Numerical Example . 

Consider the case where the system is orbited at 
the tether length is 100 Km and the loaded Centaur mass is 19,109 Kg. 
Of this mass. 5009 Kg are payload. 

For these conditions the tether mass is (See Appendix 1) . 


+ 400 Km, 


M 


T 


0.140 X 19,109 “ 2675 Kg 



PRECEDING PAGE BLANK NOT FILMED 
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Me find for this example AV - 2235 m/scc , AV^ - UA8 m/»ec 

P 

Using the Centaur data - 10870 Kg, - 3230 Kg, - 444 aec. 
the payload mass then is 




pay 


4935 Kg 


and the loaded OTV mass is 


- 18,992 Kg 

These are indeed close to the assumed values. For eomparinon, If the 

tether were not used, one trould need AV^ * 2398 m/sec, AV^ - 1456 m/sec, 

giving M “ 4356 Kg. Thus, the tether system allotjs a 13,3% increase in 
® ® pay 

payload to GEO. 


The same calculation was repeated for L = 150 Km, 19,412 Kg 
of which 5312 Kg are payload), “ 0 . 711 = 6910 Kg and, on account 

of the higher tether mass to be rewound and stored, M^p ■= 13,000 Kg. 

The results are now 


'^TOTAL “ ^23,322 Kg 
X ® 32.68 Km 
£ = 105.76 Km 

x’ = 6.56 Km 
^flN.SH "" 23.5.8 Km 


M.j,g * 23,910 Kg 
x' = 10.76 Km 
Rj, = 4,892 Kg 
jtpg* 28.51 Km 
hgH p “ 200.6 Km 


Thus, 150 Km is still feasible with a full Centaur payload, allowing a payload 
increase of 20.5% over the unassisted Centaur. However, the tether system is 
now bulky and heavy enough that rewinding operations may begin to be cumber- 


some . 
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3.5 Estimation of Pallet Maeaes for the 100 Km Caae . 

From the masses and c.g. locations of the previous examples, the 
tether tension can be Cilculated. After payload release, but before 
significant rewinding, we find 


F » 980 Nt 

whereas immediately after Shuttle detachment, 

F - 530 Nt 

A rewinding velocity of 1 m/sec is assumed. This should cause librations 
of no more than 2-4" amplitude, provided appropriate damping and terminal ten- 
sion control is exerted, and implies some 14 hrs. for each of the two rewind- 
ing phases (under Shuttle power and ovm power respectively) . 

With these data, the power required on board the lower pallet for the 
autonomous rewinding phase is 530 watt. Allowing for mechanical losses 
and some maneuvering margin, a 1 Kw power supply is adequate. This can be 
provided in a variety of ways; perhaps the most compact for this application 
would be a H 2 -O 2 fuel cell similar to those in the Shuttle itself . The mass 
and volume of cryogenic fuel needed is minimal, and the length of time when 
cryogenics must be stored on the pallet is only the duration of the rewind- 
ing phase. The mass of the 1 Kw fuel cell can be about 10 Kg, plus about 
5 Kg for reactants and tankage. 

The rewinding motor Itself must also be on the lower pallet. It must 
also be used as a generator to absorb the mechanical power generated during 
the deployment of the tether with the OTV and payload at its end. Since 
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this operation is done while the pallet is mated to the Shuttle, the 
generated power (about 5 5CW peak), can be used to supplement the Shuttle's 
own power supply, or can he radiated from a resistor bank. Allowing for 
losses, an 8 KW DC motor-generator seems adequate; at a conservative 
25 Kg/KW, this implies a mass of 200 Kg, to which we should add another 
200 Kg for gearing to the low RPM required. 

Additional mass items for the lovie' pallet include the reel drum and 
supporting structure. The volume of the fully rewound tether is about 
2.1 an aluminum drum 1. 2 m. long with a core diameter of 0.4 m and 

end plates of 1.6 m,, using 2 cm A1 . thicknesi.' has a mass of 200 Kg. A 
similar mass can be assumed for the drum supports. 

The main structure of the pallet itself. Including its Shuttle inter- 
faces, can cofiServatively be likened to the Airborne Support Equipment 
for tht lUS vehicle, which has a mass of 4160 Kg. After adding the items 
just discussed (power, mo tor -generator , reel and reel support), the lower 
pallet mass comes to 4975 Kg. Thus, even allowing for 10% growth, the 
6000 Kg used in the calculations seems conservative. Regarding the upper 
pallet, its main features may i.gain be likened to those of the ASE, except 
that, since launch loads need not be absorbed (only the approximately 1/20 
g gravity gradient force), it must be possible to lighten its structure 
considerably. Some attitude control propulsion should be added, mainly 
for control of rotation about the tether line and of out-of-plane oscil- 
lations; no estimate of these needs is available, but it is unlikely that 
the required thrusters and fuel would exceed 500 Kg. Altogether, the 
figure of 400 Kg for the upper pallet appears reasonable. 
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3 . 6 Estimated Economic Per lormance . 

It Is clear that a detailed assessraent of the econoaiea of adopting 
the scheme under dlscusalon %x>uld require a much more thorough design and 
systems study. However, some preliminary considerations can be advanced 
St this point. 

First, the Initial development and deployment of the tether system 
requires some up-front investment. Since only a medium level of technology 
is involved, an R & D and procurement cost of $40 M can be estimated. Tc 
this we roust add the initial launch cost; assixmlng the Shuttle flight can 
be shared, the 13,000 Kg tether system would displace cargo revenue of 
about $18 M. 

Let Cq be the cost per Kg for tranaportation to LEO ($1000/Kg for 

the Shuttle) and the procurement cost of OTV (estimated at $50 M 

for the Centaur). Let also , M - and M ^ be the OTV mass, pay- 

OTV pay,u pay 

load mass with no tether used and payload mass with the tether system. 


Then, the costs per Kg of payload to GEO without and with tether are 


c 


w/o 


c 

o 


a + 


^OTV 

M 

pay,o 



pay,o 


M 

* % * M 


OTV 


pay 


^OTV 

H 

pay 


) + 
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The cost saved per flight due to the extra payload allowed by the 
tether is then 

^^v/o~ ^pay " ^‘^o^^OTV ^OTV^ ^ ■" 

^ p«y.o 


and, denoting by the initial capital invescment, the nuzaber of 

flights required to pay back that investment is 


cap 


M 

' o OTV M 

pay.o 


Using - 40 + 18 - 50 M$, c^ - $1000/Kg, - 15,000 Kg, 

^OTV ^ ^ 13.31 payload increase, we find 


N - 6.9 


which indicates a very rapid payback, and justifies Ignoring discounting 
considerations at this stage. Other issues Chat need a deeper examination 
are Che possible increase in missicn^ support costs due to Che added 

complexity of the transfer maneuver, and the impact of this maneuver on 
the overall Shuttle flight costs. Some compensation may occur due to 
the reduced deorbiting AV needed after the tether release. 
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4 . P latform-Based Intfermcdlate Tether Systems 
4 • 1 Introduc t Ion 

t\s shown in Appendix 1, the mass of a tether with a given maas at 
its end Increases about quadratically with the tether length up to aome 
250 Km, after which, even with an optimally tapered cross-section, the 
mass increases much faster. The numerical examples of Sec. 3 showed that, 
for pavloads consisting of a fully loaded OTV of the Centaur or lUS type, 
a free-flying, re-windable tether that uses the Shuttle as reaction mass, 
is limited to about 150 Km in length. Beyond this length, a larger 
reaction mass is necessary, with a means of restoring its orbit after a 
launch, and rewinding becomes undesirable. In this section we consider 
systems of this type, anchored to an orbiting Space Station. Insertion 
of payloads into a LEO-GEO tran fer orbit is the mission studied in detail, 
however, other missions may be possible for a Space Station-based permanent 
tether facility, including capture and release of higher near-Earth 
satellites for inspection and repair. 

4 2 Tether-Assisted Insertion into GEO Transfer Orbit 

The system to be considered can be summarized as follows: 

a) A Low Earth Orbit space station is assumed to have a radial 
outward tether deployed as a permanent facility. It must 
also have some electric thrusting capability (over and above 
that required for drag make-up) . 

b) This tether is restricted to lengths below 300 km, in order 
to keep the tether mass from becoming dominant for its own 
tension. This length aleo provides a reasonable extrapolation 
of already planned tether technology (-100 km) . 
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c) Payloads (attached to an OTV vehicle, such as Centaur or lUS) 
are delivered by Shuttle flights to the space station, and are 
attached to a sliding "ferry" for transportation to the other 
end of the tether. The ferry must have a braking system, a 
radiator for disposing of the brake heat, controls for speed 
and some power generation capacity for return. 


d) After release from the tether end, the 0T\' engines are fired 
to supplement the velocity up to that required for insertion 
in a Hohmann ellipse leading to GEO altitude. circularization in 

GEO is made with a second OTV firing. 

I 

Let L be the tether length, M the mass of the combination space 
platforra--deployed tether and the orbital radius of the platform before 

pa\load deployment. After deployment, the payload is at a radius Rp = R^EO* 

L — j ^ where m is the mass of payload, OTV and ferry, while the plat- 


form sinks to R 


PL 


R 


LEO 

R. 


M . n. 


The velocity of the payload just 


;p P 

after release is \Ir r and after adding a perigee impulse it becomes 

I LEO LEO ^ 

the perigee velocity of the transfer ellipse, with apogee at namely 



2 R 


GEO 


* ^GEO 


Thus 



ti Rp Rqeo 


If *^150 ^LEO 


or, in dimensionless form. 


t\\ 


C, LEO 


2p 

f (f ^ 


f 


( 1 ) 
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, _U GEO . - , ^ A L 

\.,LEO = . “ \eO ' " “ ■ “ \eO ' 




The usui.1 expression for the Hohmann transfer is recovered for f « 1 , 
At the apogee, the circularization Jinpulse roust be 


or 


2R ^ P 

Av » V - !:i P . where v - -r 

» ' 'geo % ^ 'geo ^ 


AV 


" = 1 - , 


:,GEO 


V f + p 


( 2 ) 


The platform mass must be large enough to prevent too low a platform 
perigee after release; as shown in Section 3, this perigee is at 


(1 


Av . 1. 

l+v-’ 


^,Fl" 


Av. 


2 - (1 - -^) 2 
1+v 


\eo ~ 1+V^ \eo 


An example of calculations for this system is showi in Table 4.1. 
The space statior, is taken to be in a 400 km orbit (R|^gQ ” 6770 km) , 


Vhi.l= 


42200 kra. The values of v.,.„ sho^vn are those that would 

MAX 


(3) 


give a 250 km piatfoms perigee; a reduction hy 1/1.5 is assumed for safety, 
and is given as the v adopted (heavier platform) . The tether mass is 
calculatec for tapered Kevlar Aramld (p “ 1.44 g/cm^, o •= 1 .397xl0®NT/m^ , 
safety factor » 4). The payloads and inl'cial OTV loaded masses are for 
an assumed Centaur vehicle (structuraJ mass = 3230 kg, propellant mass •= 
10870 kg, exhaust ^elc Ity “ 4355 m/sec) . No orbital plane change was 
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considered. As the table shows, there is a 20% gain in payload for a 150 km 
tether, and a 38% gain for a 250 km tether. The mass of the wire itself 
varies from 0.14 to 1.27 of the maximum end mass. Since this mass is 
of t.. rder of 25 ton in this example, the maximum tether mass (for 
250 km) is about 32 tor. The platform mass varies from 5.3 to 14.7 times 
the end mass (as a minimum); i.e., from 130 to 370 ton; presumably, this 
would include the empty Shuttle attached to it. All these figures are 
reasonable, and appear to be within the scale of the contemplated Space 
Operations Center, or expansions of it. 


Tether length (km) 

0 

100 

150 

200 

250 


0 

0.01477 

0.02216 

0.02954 

0 . 03693 

MAX 

0 

0.273 

0.167 

0.120 

0.0937 

V (adopted) 

0 

0.182 

0.111 

0.080 

0.0625 

M(Tether) /m 

0 

0.140 

0.356 

0.711 

1.269 

M(Platf. only)/in 


5.35 

8.65 

11.79 

14.73 

f 

1 

1.012S 

1.0199 

1.02735 

1.03476 

AVp (m/sec) 

2398 

2233 

2133 

2035 

1938 

AV (m/sec) 

1456 

■ 

1447 

1442 

1437 

1432 

(m/sec) 

3854 

3680 

3575 

3472 

3370 

\ (kg) 

4409 

5009 

5312 

- 

5684 

6076 

Mq r loaded OTV) (kg) 

18,509 

19,109 

19,412 

19,784 

20,176 


) Centaur 
OTV 


Table 4.1 Perforaance of tether-assisted LEO-GEO system. 

R, „ * 400 + 6370 km„ 

LEO 
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4 . 3 The Platform Propulsion System . 

The climbout of the ferry would lower the platform c.g., and the 

release of the payload/OTV would send the platform Into an elliptic orbit 

with perigee well above the atmosphere. A propulsion system is required 

on board the platform to restore its orbit before the next launch. The 

thrust can be applied either after or during the ferry excursion. 

‘.ercury bombardment ion engines have been developed to the point where 

confident performance and mass estimations can be made. Byers (Ref .4.1) 

presented a methodology based on extrapolations from existing thrusters 

which can serve as the basis for our analysis. Specific impulses (I ) 

sp 

from belnv.7 2000 sec to over 4000 sec are possible by adjustment of 

voltages. Very low values of I lead to high propellant resupply rates, 

sp 

as well as to low efficiency of the thrusters. On the other hand, very 

high I implies high power requirements, with attendant mass increases, 
sp 

We present next a study to determine the appropriate specific impulse for 


our application. 
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Thfc input power to a battery of ion engines operating at exhaust 
.’elocity c, with propulsive efficiency and thrust F is 


P = Fc/2’- 


(4j 


and in terras of the velocity increment IV to be imparted to a mass 

'1 in a time t, , 
b 


P = Mci.V/2’- t^ 
P b 


(5) 


In Appendix .1, an expression (Eq. (28) of that Appendix) is 
derived for the IV required to re-establish the orbit of a space 
station at R, „„ after release of a payload from the end of a tether 

L h J 

line of length L: 

IV = 2.352 V )(|' ) (6i 

^EO '‘Total 

where m' is the mass released (OTV + payload), L is the tether length, 
and v^ is the circular velocity in LEO. If the engines operate after 
payload release, the mass to be accelerated is M = • Also, 

Eq . (5) gives the average power during orbit recovery, but, as shown 
in Appendix 3, the thrust must be applied in a modulated fashion. 


F - Fq “ "I 


(7) 


where 6 is orbital azimuth from perigee. This leads to a ratio 


^ K/', 

max ' max 5/2 

“ <jF|> 1.176 


Therefore the peak power required is 


P 

max 


f ^ - a - ? 

2 ^ ^^TOT 


) (- ) — 


(9) 
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Some modification is needed if platforji thrust is also applied 

f 

during the ferry cllrabout phase, but since rr is tjrpicaliy - 0.1, 

TOT 

the impact on P is minimal. Notice the small sensitivity of P 

Tn3X 

to H , and the proportionality with tether length, 

platform 


The .-’mount of propellant (H ) used follows from Eq. (6) : 

£ 


^ “ 2.352 m'(l - f' ) (|- ) ^ 


( 10 ) 


The propulsive efficiency of ion engines increases as the 
specific impulse at vrliich they operate is increased. 

In general one can write 


'cD 


2eV 


!-*■(: 


LOSJ 


tn^ c 


( 11 ) 


where n is the power conditioning and distribution efficiency, 
cu 

V^oss thruster power loss per ampere beam current and e and 

arc the ion charge and mass respectively. From the detailed 
analysis of Byers (Ref .4.1) one can use for existing and near term 
mercury ion engines at 0.95 propellant utilization fraction the 
values 


n = 0.752 


''loss' 1” ''“1“ 


this gives 


0.752 

L “ 

P 1 + 1.282xlOVc^ 


( 12 ) 


(an almost equally good fit can be obtained for the more physical 
value -= 150 Volts if n is raised to 0.765). 
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The cost per niysion Includes some compom'ints that arc sensitive 

to the choice of exit velocity c for the ion engines. These are 

(p) A recuriinc cost Cj^j , where is the cost of 

E B 

mercury per Kc (in orbit) 

(b) Kon--recu’'riPS costs, mainly the cost CM of the powo- 

P P 
s ^s 

system, v.’hcrc M = a P and a is the specific mass of 
^ s 

tV»e power system (Kg/watt) . Other non-recurring costs Chat 
may depend on c are those associated with the ion engine 
hardware; higher specific impulse implies smaller fuel 
tanks and other fuel-related components, but larger power 
cenditioning and power-related components. Overall, Ref. A.l 
concludes that the engine system mass is insensitive to 
specific impulse, so we omit this from our discussion. 

V/e are thus led to choose the engine specific impulse c/g by 
minimizing the partial cost 





g 



a P 

max 


(13) 


where N is the number of reuses of the power system. I’sing Eqs. (9; , 
(lU) and (11) this can be rewritten as 


<P 


aC 


1.063 




N n. 


(c 


+ — — — — 


(lA) 


where 


2 \oss 

'£ m. 


1 T1 

1 w 


1.063 c ^ 

ps 


(15) 


o 


( 16 ) 
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Differentiation gives 


f> = '\ f c ^ 4" c ^ 
OPT \' )L a 


(1 


i.e., ('qpj. ■■ Cj, for no reusability, but c^p ~ N for many reuses, 


An estiniatc of the cost of a multi- hundred Kv; solar array can be 

obtained from Ref .4 . 2 , where detailed design and costing is performed 

for several types of arrays in the 400-600 Kw power range. The 

lowest cost (for lov; concentration ratio GaAs cells) was found to be 

326 $/Avg,Watt, of which about 90 $/Uatt corresponds to launch costs. 

The array specific mass v:as also found to be about 10 Kg/Kw(B0L) . For 

an assumed ratio of average to BOL power of 0.85, this leads to 

27,70( S/Kg array cost (7630 5/Kg for launch). 

The high array launch cost just me’^tioned is related to the special 

arrangements for pressurized Shuttle bay stowage and self-deployment. 

By comparison, supply of mercury propellant to the space station is 

likely to be a simple operation: we assume a cost of mercery in orbit 

of Cjj = 2000 $/Kw (including a comparatively minor allowance for 
8 

purchase price) . 

The power system specific mass cx includes not only the array itself , 
but also other componentc . such as gimbals, regulators and battery system 
for eclipses. For the first two items we follow Ref .4.3 and assume the 
following masses: 

Gimballing system 4 Kg/Kw 
Regulators 5 Kg/Kw 
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4 3 


For the ba r t c r J o;; , wo aEsuriu? type, with energy density 

J7 ' att hr/Kc ^nd clia ri^e-dlschargc efficiency of 0.77 (Ref. 4.3'). 

The lota] energy storage aoeded can be calculated using the thrust 
profile of be;. '!) if tlic eclipse time is specified. Tlic worst case 
for shadovv’ing occurs vdicn the sun lies in the orbital plane, and 
gi\'CS a slindov.’ tiir.c 


b.' A/ 


' n 

LEO . -I, T. , 
; Sin (- ) 


\eo 


(18') 


where is the Eartli radius. 

An additional consideration to be made pertains to the relative 
location of the eclipse zone and the orbital perigee; this is 
important, since, according to Eq. r7>, the perigee power demand is 
only 1/5 of the peak demand (at apogc c) . Apogee for the perturbed 
platf orm orbit occurs at the location of payload relea.se, and one can 
in principle place it at orbital noon to minimize energy storage. 

Since the fuel cost of moving the payload within a GEO orbit is small, 
it seems reasonable to assume such a release strategy. With this 
assumption, the mean power demand during eclipse is gix'^en by 


<P .> 
sh 


max 


sh/2 


sh 


where 6 , = 2 sin 
sh 


0 

-1 


2-3 cos e| , 2 0.2216 - 0.6 

5 ^ 


( 19 ) 


sin 


(— ) 
\eo 


R^/R^Eo orbital arc in shadow. A similar 

calculation can be made for the capacity excess during the sunlit phase 
of an orbit for an array dimensioned for peak power. 
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Some results of these calculations are shown in Table 4.2. 


Table 4.2 


Orbii u]tiLudc (i'Ji) 

100 

200 

300 

400 

500 

60C 

Max. Uliaci '.’ timo (r'.inutes) 

38.37 

37.26 

36.56 

36.08 

35.73 

35 .46 

<P , >/P 

Ell iiiax 

0.1353 

0.1279 

0.1232 

0,1201 

0.1181 

0.1166 

Idtnl Btorngc rc^ (Kvh/K-.A 
max 

0.08652 

0.07943 

0.07507 

0.07222 

L .7033 

•' t6891 

Excess canacitv during sunlight ^ 
storage required 

3.009 

3.516 

3.976 

4.402 

4.799 


The last row of Table 4,2 shows that no extra array area is required for 

battery charging. The storage required is only weakly dependent on orbital 
altitude. Using the value 0.07222 (for 400 Km orbits), and including a 

battery efficiency of 0.77, the battery mass needed is 

B attery mass ^ 0 . 07222x1000 (Watt hr/Kw) _ . ^ 

Peak power 0.77x17 (Watt hr/Kg) " Kg/Kw 

Thus, including the array, gimballing, regulators and battery system, 
we arrive at a power source specific mass 

a = Q 85^^ ^ ^ 27.0 Kg/ (Kw to engines) 


where the factor of C.85 accounts for the extra "BOL array area required to 
accommodate cell degradation. 
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For the case where the Cec.ff’ur OTV is used, the mass released Is of 
tlie order of 20,000 Kg (See Sec. 3), varying slightly with LEO altitude 
and tether length. For power estimation purposes, the ratio will 

be assumed to he 0.06; this is compatible with a safe platform perigee 
height, and in any case, is an insensitive parameter (Eqs. (9), (10)). 
Finally, we choose a total ion engine firing time of 14 days; as we will 
see below, this is about twice the ferry roundtrip time adopted in this 
study, and should therefore set the maximum mission frequency for the 
tether system. 

With thes-? parameters, Eqs. (17), (9) and (10) read 


c 


on" 


V 1.282x10® + 2.289x10*^ N 


,>(R0L) 3.8S4j;]0 

^lAX " "OT" 





8 


= 44200 ^ 

\eo 


V 

c 


c 


Table 4.3 siiows calculated results for a 250 Kra tether 


N 

1 

10 

30 

100 

300 

— 

‘^on 

(m/ sec) 

11,420 

12,290 

14,030 

18,900 

28,550 

”sp 

^ OPT ^ ^ 

1 , 166 

1,254 

1,432 

1,928 

2,913 

■ 


0.380 

0.407 

0.455 

0.5533 

0.6498 

B 

P (KW) 

max 

331 

332 .3_ 

338.8 

315.5 

483 _ 

IB8W 

— ^9 

390.9 

398.6 

44i78 

568.2 


(Kg/mission) 

1,095 

1,018 

891.9 

662.2 

438.4 


Table 4.3 Optimized propulsion system parameters as a function of 
number of reuses, for a 250 Kjii tether and a 400 Km LEO 
orbit . 
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and a 400 Kn platform orbit. These results are ineensltlve to orbital 

lieight, while and scale in proportion to tether length- As 

g 

shown, the range of specific impulses from 1500 to 3000 sec is optimum, 

depending on reusability. The tradeoff between power and propellant 

mass is apparent from the last two rows of Table 4.3. 

If we adopt 1 “ 2000 sec (optimum for about 115 reuses), Byers' 

sp 

analysis (Ref. 4.1) can be rather directly applied. The accelerating 
voltage and net voltage (including the decel electrode) are 2000 and 443 
Volts respectively. The individual thruster diameter was selected at 
50 cm. The results are summarized in Table 4.4 (for a 250 Km tether). 

Table 4.4 Platform propulsion system characteristics 

Type 

Diameter per thruster 
Specific impulse 
Thrusting time 
Thrust per unit 

Thrust power per unit (including 
distribution losses) 

No. of thrusters required 

Mercury mass per mission 

Solar array power (EOL/BOL) 

Thrust system mass (thrusters, thermal 
control, power supplies, interface 
module structure, etc.) 

Solar array mass 

colar array gimballing mass 

Solur array regulators mass 

Battery system mass 

Total propulsion related mass 


Hg ion bombardment 
50 cm 
2000 sec 
14 days 
0.546 Nt 

9.56 Kw 

41 (+ 4 extras) 

632 Kg 

384.2/452 Kw 

5464 Kg 
4520 Kg 
1808 Kg 
2260 Kg 
2113 Kg 
16,165 Kg 
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4.4 The V r i r y Di l ve Sy.s tfr:'. 

Ill this SLCtion i.v caJctilal,e Uic required power generation and pov’cr 
dir,£.j ; ..t ion (..parities of the ferry vehicle tliat transports the payload 
and O'! V to the end of the tether line. 

Let i.i he the outbound travcllinp, iiiass, made up of the OTV , the pey- 
ioad and the returnable ferry 


m*=Ei'-4-Mp ^ (21) 

vdiere n ' ’’‘l OAU mass. L'hen the ferry is at a 

dir.Lancc y from the lower platform, its distance from the (moving) 

overa)] center of mass is y-y^ “ (1 - v) y - v L , where 

eg 1 


M + tn 4- M 
P 1 


M„ M 

T p 

N +m + M,^ * M +m + M_ 

pi p T 


( 22 ) 


and M 


Hence , 
dv 

3 l-> 

dt 


and 


tlie 


is 


M,j. are the lower and upper platform masses respectively, 
raochanical power being generated when the upward velocity 


^up (23) 

Notice P <0 when y < r — L ; i.e., external power must be 

up 1-v ^ 

supplied to reach this point (at which time the ferry is at the overall 
c.g.). For y > — L, power is being generated. 

A similar expression applies for the return trip, when the travel- 


ling mass is the ferry alone: 
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platform. 

(c) Radiation of all the generated energy. 

Regarding option (a), we notice that it would allow elimination of 
a separate pov;cr source for the return trip, such as a solar panel 
(supplemented by batteries for eclipse times). Thus, the option can be 
assessed by comparing the required battery mass needed to that of the 
displaced power supply. 

For the power supply, if one is used, the BOL array power required, 
assuming 85^ degradation at EOL, 75% DC motor efficiency, 77% battery 
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efl.ciency and 36.1 min. shadow time/36.4 min. sun time (400 Km orbit), 

must be 

p . !downj^ J_ . 36,._1_ V . 2.68 P 

sa.BOL 0.7 5 ''0.85 0.77x56.4'' down, max 

The mass of this power system is then 

M . P + 

ps °*sa sa,B0L 3 

where, following Ref. 4.3, a *12 Kg/Kw (blanket) + 5 Kg/Kw (regulators) 

Sfi 

+ 4 Kg/Kw (glraballing) - 21 Kg/Kw. 

Also, E is the energy to be stored in the batteries for the 

Do L L 

eclipse time, and B ■ 17 Watt h/Kg is the energy density of the assumed 
batteries, so that 

^att _ ^dovm,max 36 .1 1 ^ ^ „ 

e 0.75x0.77 60 0.017 down, max 

Altogether, then, 


^8 •= 118 Kg/Kw 

^down,max 

The mechanical energy needed for ferry return can be calculated by 
integration of Eq . (24). This leads to a battery mass estimate of 


- V^) 

^batt 0.77x0.75 6 


and comparing to Eq . (24), 


M 

batt L 1 

* V 0.77x0.75 6 

down, max 


1-v’ 


- V, 


l-v'-V.^ 


( 25 ) 


(26) 


(27) 


(28) 


(29) 


(3) 
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For a constant velocity return, L/v is the return time. Using 
B 17 Watth/Kg, * 0.02, = 0.01, we obtain 


M 

* 1.20 t(days) (Kg/Kw) 

down ,max 


(31) 


Comparison of (28) and (31) shows clearly that, except for 
unreasonably fast returns, storage of power requires much more mass 
than direct generation via an on-board solar array system. 

Regarding option (b) (propulsive use oJ power generated) , 
a simple calculation will show that the contributioi\ to the required 
AV for recovery of orbita] platform is too small to be worth considering. 
Tlie thrust that can be generated with a power P is F « 2r|^P/r. . Also, 

e 1— 

M 

TOT 



Thus, using for the power P - 0.75 using Eq. (24), 

we obtain 


AV 


2.Z5 


V 




1 

(-T + V 
m 2 


V ~ 


K) 


( 32 ) 


For values of the variables comparable to those used in other parts 
of this report, this AV amounts to less than 1 m/sec. For comparison, 
typical required AV values for platform orbit recovery amount to 50-100 in/ sec. 
Therefore, it does not seem advisable to include electric thrusters in 
the ferry for primary propulsion. On the other hand, one can expect a need 
for attitude control and out-of-plane libration control of the ferry; these 
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nocHtbi have not been quantified yet, but the ready availability of the brake 
electric power may make it attractive to parfonn theae taaka with ion 

thrusters . 

Following the above arguments, only option (c) (radiation of all 
(or most) of the brake power) remains. This would appear to pone no 
special problems, since the power 1 b in electrical form and c^n be radiated 
from resistive loads whole design temperature can be quite high. 

Consistent with this design concept, a aolar r.rray power supply is 
tuHided for the return trip, with a mass/power ratio of 113 Kg/Kw 
(Kq. (28) ). The peak power needed depends directly on the mass of the 
ferry, whiclt has not yet been determined precAaely. For an OTV-payload 
combination of 20,000 Kg mass, a preliminary eatlmate la • 30C0 Kg 
(DC motor-generator, controls, guiding pulleys, OTV attachments, trusalng) . 
Following Eq . (24), for a tether length of 250 Km, and a ferry apeed of 
1 m/sec (return time - 2.89 days), this gives a paak mechanical powar 
requirement of 3.b Kw, and therefore a aolar array-battery ey»t«s: mass 
of 670 Kg. This can be easily reduced, however, by operating the ferry 
at a lower speed near the end of the tether, where the gravity gradient 
force is largest. 

4 . 5 Dynamics of the tether system during ferry transfer . 

the tether system is permanently deployed and payloads, with 
their orbital transfer vehicles, have to travel along the cable, new 
dynamic effects may arise which have not bean dealt with in the literature. 
For instance, the ascent velocity v of a ferry of maaa m gives rise to a 
backward Coriolis force 2mnv, which leads to oaclllatory in-plana motion 




of both, the ascending maas and the two end masses (the main platform 
below and a small terminal platform above) . At least some of the modes 
of oscillation have the feature of a rapidly increasing frequency as the 
distance between two of the masses approaches zero; there is therefore 
the potential for a wrap-around type of Instability when the ferry 
approaches the end of the tether. Similar effects can arise due to the 
tether elastic Vty, and here the risk Is that of greatly enhanced tension 
due to dynamic ''ffects. 

Two lines of attack !iave been followed in this problem. On the one 
hand, an analytical theory with some simplifications was worked out, first 
for the in plane oscillations (Appendix 4), and then for the stretch 
0 c i list ions (Ref. 4.4). This separation of the problem into two individual 
problems is allowable because, due to the linearization used, the tvro types 
of motion decouple to the order retained. Independently, the same problem 
was treat -d at the Smithsonian Astrophysical Observatory (SAO) using 
numerical methods which bypass the need for many of the approximations used 
in the analysis; these SAO calculations will be reported separately. 

Both studies reached very similar conclusions : for a tip mass of the 

order of 10* of the platform mass, enough tension exists in the tether to 
prevent instabilities and maintain oscillations within fairly small bounds 
for climbing speeds of the order of 1 m/sec. The only time when a divergence 
may occur is the terminal approach phase; and even there, careful speed 
control in that phase can ensure a smooth maneuver. The detailed analysis 
for the in-plane case is given in Appendix 4; a WKB approximation was used, 
together with an inner-outer matching process near each end of the climb. 
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Thc; ttubKcrijit 2 hat. )xcn \o indicato condU;jlunt» 

ciftior cp, 'Hc-otion oi . FoJ th*,j iiiU!0?+. oirl'ii t (bt'fux'o » t)K> 
aprvyctc. vt>:Ui»”.r1 y J s; 


Q 

V(cipc»y:'». j^) >- V (iiprajf' j^l g — “ ^'g 

GS 


(f>) 


iind Fq. (.v) cau be wc'cHficd to calci\j£itv the pc"igfie R_ of thit 


aacoat; orbit: 


2 

.1 (■ ~~ \ ■.- IV — ' « Av ) 

‘ ‘V,- \ “cs ^ 


a) 


Cnfo r. ir; i-.o d<tcr;,-:nc;1, the velocity at perigee can be esiv'enrod 
3 

(troKi <;oi)x.crvitij.on of a;iga7.r.i' nfioii'cutinn) as 


V « V(fipo.,) « (V Av ) 


( 8 ) 


TJiisi vc\lor,ity coni e, ins, in cjenora,l,a propu .Ik. ion -derived incixment 
AV„ , applied at or luvnediatcly after rc.lcaae. ^The veloc.ity of th<’' 
end of the low tether is therefore 


V . » V “ Av 

tether P P 

end 


( 9 ) 
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.Tvi U’.'o ojbit-iJ v< lo'^.lV.y V^. of thr at. 

(or, »i* ro ju- oo.lt-c.lv, ol the- o.t »-,j lua.! center of t\.c tolhar i»ayJo-fcl system) 


ir. t1'»' rc.Co; <’ 


;,),n “ / ~r- ” ‘'V, ■ k' 

’u: ^ 1 


( 10 ) 


fro,.! vhir)i r,,._ Cit:' »k; oaXct.vly.tcn cc-sUy. FinoHy, the 3ow-13tirth l ethe r 
Icsol-h if- 


Pig. 5.2 moA Table 5.1 tliov oalciilated r««ilte for tbo com of 
LV ^ « AVq ■ 0. If M i^poM tbo roquiroMHit chott in com of docldng 
foiliiro» the poplood ond the Xowor plotfoca of •’^90 tothoar should 
coBdossottS ogoio oftor ou iutogor nwabar of oi^i'Co» chon tho period Pa 

(T In Fig. 5.2) anst bo o rotioiiol froction m/n of o dop (m* n iacogoro) . 

Tfano* approprioto voluos of Pa for loo a and a, ore 

1/3 day - « hr,, 3/8 day - 9 hr., 2/5 day • 9.8 hr., 1/2 day - 12 hr. 

As shoim in Pig* 5.2, a period of 1/3 day inpllM an upper tathar 
length of over 10,000 fti, and a loaar tathar length of about 1200 in 
fcoB a lev larch orbit at 1200 la m vaU. lacraasing the period to 
1/2 day lovers the length of the nppar tathar tc about 8000 la, but it 
alM raquiTM the lov Barth orbit to ba at sene 9000 la altitude, vith 


a 1800 la tatliar. 
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Or)>lt i!rom LmO to UHO 



PI6DSE 5.2. 


TWO-TIJTIIKR SYS’jnM ClLMlACTKlUS’fJCC WITH AVp - « 0. 




Sf 


Of POOR Q»**I-ITY 
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TiMi VtlI7 iMIt tSCiMHM fOT thiS MtUMllf fUlt* 

h t w y, mi 4 will iMvt to ko p ■—■■■■ tly ^lopoi. Thovo io • otioiit ineon- 
tl»o to OKploit tko oBpoMittol difowiffineo of moo oo looftli Op opplplag 
port&ol Mpololoo thmiot ot tlio oaio of tto trooofor. 

Tho totlior MOO colcolotiM io cMtolood io Appoodix 1. It io of 
iot^nrast thot tho Moh ioofor CIO totlior io of oboiit tho oom ooob. 

Yho offoet of iatrodoeiof both porijoo ood opogM firiopo (4V^ ond 

voopMtivolp) MO aoot loMOtifotod* Tho rooolto for o vido roogo of 
poroMtoro oro liotod io Tohioo S.2 tdnmi|h S.5. for tho ooom of tho 1/3 
ood 1/2 dop poriod* tho rooulto oro oloo di^lopod tf^ib^Leollp la Fipo. 5.3 
ood 5«4. Tho offoeto oro poaorollp m folloM; 

(a) incica; i 4/<j at, i;oiutt ium. incjrofiOOF iiio altitiuf''*-* oi j#i:C, 

and t-Vic* ;iovA*r totla'i JcnqUi, h. 

(b) Incrcasiny AV^ at c€ip(;tpnt the altittnU; of J >:o, 

roiT low Av , ancj t)UP'.i!i of AV rraul t in a shoJTl or Iova r ti'tli. r, 

• VI 

but Uir rvVt.ifU! in trnu at hifii) vaXucii of AVp WO w/ficjf) , 

(c) hi'. d.j Bt ui cvf) heforo, ihf‘ Icnqib H of the upi> )• tetl-ujy ir: i*nr f fw c- 

tod by cithev A\^ or but i.*? red ;ooc1 .if i-i>e p(jr.ioA ,ir. SiJlOi.'v-d 

to ,lncvi:ai!c. 

<d) Tor each t ranrfer I into and each valite of Av_ , there fn n 

r 

w«in fet v'lich t>»o lower Karth orJ:>it boco!r.a« too low (a liiuir 
of 200 Kr: ar.sv’nyd ncjo). Sindlarly, for each there ic; 

R m:i.nintu*n AVp for the uat«' reouen. 

(e) Tiiu loTifjth of the 3r*.V' r tether can ho rertticed to zero hy inereoe.i 
AVp fc»r fcel AV^^ . Ths effect c>f AV^ on h is minor. 


fMilftfMAL f9M It 




Ku f'£ it iJiniibinotift.i wUich c»o«ia l*c* u^ofyJ , w© ■<»« in 

4 iliat# #r©r.i »• Km m ui»irf;i it Jov<r itttoer of Icniith 390 1*» *'.n0 

iiupp?yi»0 H vcJc city iiiCKeiatmi' /^V^, ,1600 aCtt © KoloaitOf « payload 

o©.n be. put i’.iit' tt tifu»u»fe» f’l.).ipn' to ejitptui'© by tl»o lo»ir©r jmd ol 

• CKO U'tl.iU* of K»'» .length# .if ca apo^uu volooity iucresuoni 
AV^ « VI?*:. h./noc it* ip >l.ied pi*ior to dowKiiig. If tlooklng miothor 

cati bo ii*»uk attoi. ont? <l«y (tv?o ortit?. of tlio }3ayl.oad). Nolioo 
that fot thio Iwoi.* totlior length# ito r»qs« oan l*«* of the oird«i* of fcUt’ 
payload tivma. 







































o 

r* 


fl* 

Jl 

4» 

T1 

v4 

<•) 


»4 rf! 
It H 


k 


o 

14 


i>.4 

0 

0 

1 

4 -' 

f‘4 

<! 


V*' 

js 


a» 

« 

w 

•ri 

M 

4J 

c 

14 


















M 


ORKSiNAL PAGE IS 
OF POOR QUALITY 


Orbiif.'* |>r.y1u rhot ’ lr>un of ihr i l OolgonB 


Jn the chlcul»itlt>nt5 r.o fair \m have implieit-ly very ha*Avy 

ijlafcJ'oi.rt'O both l.j» I.no out! in G!‘0. If Um raann of tlio 7,%0 platform in 
<l(D*n^ tintfO }»vt that; f»f the* .fihultla orbitcr (dj*cUo<1 to a light frao -flying 


teUi'^J- facility, tiio £iilio Vfy-^C).ul 

(11:1 for a .'t.C 'i’on’ia puylo; '■;) • I’ho rq.wlt nuy be in cxcfr.f;iivo lowering 
of. the jyc)»t* rolc?’*?'^^ perigee. 1« thin ocicfcion we consider this 


effect* while r.lill ftSMurtiiug a inseeive GEO pletforfii. 

I'he J»,ev 7 gecji.otricaj, firrangomeat fc^r the lower tether is shown in 
flu, 5 ^ 5 , the orbital confre;.* is at # given (l»f .5.1) 



h 


Flgitvs 5.5. Geometry for a finite lower platform mass. 


l(rf.5.ll Study of’oaria’in Twobiiiqties Using hemg Ch'hitlng Tethertj 

by Giuseppe ColomTx*, B'iital Report on grunt MIM3-8008, from the 
»AO to lipSA. Feb. 1901. 



OIIIGWAU 
OF ipS&B Qu/^nrY 


R- 


r 1/3 


£ «j/r| 


(17) 


aad !• loc«t«d at m dUtam^ h* * R«, ** K eh« eraimfar orbit 

h ^ 

ftorlgoa. flHia* h' roplaooi h juid R^ raplaeoa in otir proviouo 
•lUilyala (Iq«. (10) , (11) ) . Tbo nav Boiat ba obtaii^ froo tha 
axpliclt lEonn of Eq. (17) ; for aacaaple, accounting only for ttio and 
naanaa Hi, Hs (Fig. 5.5), ae have 




1/3 


Mi/R* 


(18) 


M 


XiBO 


which can be aolved for . 

The parigae of the post‘*>ralaaae platform orbit can be calculated 
from Iq. (6) of Ref. 5.1, which for our case reads 

82 ■ - *tEO + ) (£9) 

The affaet of this modification is to require a longer lower tether 
and to mal» high values unfeasible (negative perigee). M an eamuple. 
Tables 5.6 and 5.7 show a ccmparison (for 1/3 day period) of two cases, 
one with a massive LEO platform (Hi ■ 5000 Tonne for H 2 ■■ 10 Tonne). 

In the first case, where only a slight perturbation is introduced to 
the orbit, a teth^ length h ■ 998 Km can be used from a 521 Re orbit, 
which becomes a 521/SU orbit after r^ease. Velocity incrmsents AVp«* 300 
' m/sec. AVq ■■ 100 n/sec are required. In the case with the light platfOTm, 
the AVq 8* 100 is not allowable, and so, for AV^** 300 m/ sec, only AVq« 0 
is possible. The result is a longer tether (1155 Kn) and a higher orbit 
(1291/858). 
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TAMLE S.6. i>i,ATFO»l IK l,EO 

Ml w 5,C0f),0>% Uj' (iilfitConn) 

M2 10»000 lif’ (ttfitclHlci) 

1* « 1/3 day II «« 10*390 Uw 



lincriee arc h in ktn 


per.! g;c v> 

nlt;I,tvide 

lU'.^ativc 


apoj^Gc altitude 

(Urn) 


perigee altitude (luii) 
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T«thar Prop^rei— and Tmtlmt 


1. Calculation of T«th«g 

For a radially daployad constant etraos tathar (atraaa ■ a, danalty " p) 
tha eroaa aaetikm A(r) ouat ba maximuiD at tha orbital eantar (approxiJiiataly 
tha c.g.) of tha orbitiag aaaattbly. Lat ba thia naximm aaetion and 

tha orbital radius out to tha c.g. Wa can than aaaily find from statiea 

that 


M.P 

A<r) . «pt— 


/ 3 . r^ 

^®LBO 



a^/ 


Cspanding and retaining only quadratic tarma (or» altamativaly, start- 
ing from a constant gravity gradient appr^ination) , 

r-I 


q r**K| 2 


< 2 ) 


At tha upper and lower end of the tether, the respective concentrated 
masaea and nuat be in force equilibrium between tether tension 
and gravity gradient force: 




3 ^max -PI- 2 < 


^“*«o 2 

] 


(3) 





( 4 ) 


where L is the stretched tether length and x la the dletance from the 
lower platform to the overall c.g. of the system. 
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EqiwtioM (3) and (4) can ba aolvad Cor and Unforeuoataly, 
axcapt for Uniting caaaa, thia aolutlon cannot ba obtained in eloaad fora. 
To faeilitata diaeuaaion, let 


Y 



3 V 
2 0 


\eo 


(5) 


V 


!!tiR . 

“PL ’ 



(6a, b) 


Than, by division of (3) by (4), and after aleqpllfieatlon, one obtains 
an aquation for ^ (c.g. position): 



1) . ,-r^(x-25) 


For abort tathare (y « .1), this has the approxlnate solution 


(7) 



( 8 ) 


For other conditions. Table 1 lists values of ^ obtained from Bq. (7): 


I 

1 

v«0 

0.05 

0.1 

0.15 

0.2 

1 

o 

0 

0.04761 

0.09091 

0.13043 

0.16667 

1 

0.5 

0 

0.05869 

0.10843 

• 0.15150 

0.18937 

1 

0 

0.10011 

0.16380 

0.21097 

0.24823 

1.5 

0 

i 

0.17652 

0.24201 

0.28396 

0.31499 

2 

! 0 

1 . 

0.25775 

0.31139 

0.34371 

0.36697 

3 

! 0 

1 

0.37446 

0.39562 

0.41393 

0.42694 

4 

i_L_ 

0.41687 

0.43608 

0.44732 

0.45531 


L 


Table 1. Values of C 
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Ttui naas of cha caebar can ba obtalaad by Ineagraelcm of l^i* (2) : 
*■■*<!» 

Hr * P Aj... /„ ««pl-Y*(ft*ldy * P A... ^«rf(YC)* «rf(Y(l-g))l {#) 
eg 


whara y " * %»a • Tha va3.ua of A la obeaiuad from Bq. (3). 

IpUv nax 

Aftar aooa coarrangaaaiitf wa obtain 


^ - 7T [arf(YO + «f(Ya-C))l (9) 

“tip 

For anall y (abort tatbara)* Eq. (8) can be uaad approxiaiatlay for Ct 
with tha raetilt 


M 


tip 


JL. 

14^ 






( 10 ) 


a 2 


jd 

1+V 


(-^) 

'W 


( 11 ) 


The laat £orm» valid roughly whan < 0.3 , indicataa a 
quadratic dapandenca of maoa on length for short tathera» where tha tip 
naaa doninatae cl^rly over the tether naaa. Whan y approachaa unity* this 
changes to a nueh stronger exponential dependence* as the mass of the tether 
Itaelf becooes doninuiit in deteminlng its cross-section. 
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T«bl« 2 Use* valuM of for noro goiitrol eomlitiono (fron 


(9). 

uaing Eq. 

(7) for O 

e 




v«0 

0.05 

0.1 

0.15 

0.2 


0 

0 

0 

0 

0 

BKl 




BEEmi 

0.4^/3 


(0.5923) 


(0.5259) 


(0.4315) 

1 

4.0002 

3.2^00 

2.8218 

2.5485 

2.3136 


(4.0602) 

(4.3256) 

(3.6496) 



1.5 

24.370 

12.195 

9.4296 

— 

7.0652 

2 

192.640 

35.713 

25.530 

20.920 

1220111 

3 

43.090 

239.68 

163.17 

130.85 

112.10 

4 

6.30x10^ 

1888.5 

Hi 

1032.2 

884.7 


T.U* 2. VUuM ot 


Tho flgisros in poron'ihoolB In Toblo 2 Aro calcularod according to 
Eq. (10), for eoflq>arlson. Those raoulC 0 are preaanead graphically in 
Fig. 1 (for Y < 1.1) and Fig. 2 (for higher y). 

for purpoaaa of ealibraticm, let ua aeatime the following propertiea 
(appropriate for Kevlar tethere): 

p ■ 1.46 g/cm* ■ 1440 Kg/m^ 

0 - 1/4 140 Kg/inn^ - 1/4 1.4xlO®N/in^ 

and alao Bg + 400 Km • 6.77xl0^m. Ue than calculate 
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3 



Prom Pig. 1 m ••• nou that tha tathar n»aa Inccaaaaa lika for 
L 200 - 300 Kn; for longar cathara In LEO, tha naaa aaealataa rapidly, 
as ahovn in Fig. 2. For L « 1200 Km, aa raqulrad for no-propulalon 
tranafar to GEO, va aaa that 200 . vhila if thin langth la raduced 

to 600 Km by uaa of partial propulalon, than 15. A eimilar 

raduction oeeura if tha working atrangth could be doubled (aaa dlscuaaion 
below) . 

If ie replaced by • 42200 Km, than 

V ~ k(^) 

’CEO " 5500 

showing that much longar tethers can be deployed in GEO orbits. 

Properties of Tether Materials . 

The single moat Important property of a desirable material for our 
application is a high specific stress (o/p). Fig. 3 compares the a/p 
data of many high-strength materials, including steel, fiberglass, boron 
and graphite fibers and the fibers known under the trade muse of Kevlar 29 
and Kevlar 49 (Dupont). The latter are clearly the best candidates, unless 
high modulus is important to minimize stretch (in which case boron or 
graphite fibers are superior). A similar comparison, this time in terms of 
the direct stress-strain curves for several fibers, is shown in Fig. 4. 
Values of a up to 3.6xlO'^N/m^ are shown for Kevlar in die form of Impregnated 
strands (360 Kg/mm^). 

Physically, these strands are made if a bundle of very thin fibers 
(diameter -12 vm), and the values quoted refer to tests made on samples of 
a few inches in length. Clearly, the probability of a flaw increases with 


SPECIFIC TENSILE STRENGTH, 10* IN. 
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H-6 


FIGURE 3 . SPECfFtC^ TENSILE STRENG TH 


SPECIFIC TENSILE MODULUS 
OF REINFORCING FIBERS 

TENSILE STRENGTH OR MODULUS OIVtOED BY DENSITY 
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ilSSSO, STRESS -STRAIN BEHAVIOR OF 
REINFORCING FIBERS 



RESIN IMPRESNATEO 
f STRANDS 
CAtTII 0*>2S4a) 


60R0N 

FIBER 


ROVING OF 
KEVLAR® 49 
ARAMtO 


TYPE HT 

GRAPHITE 

YARN 


DRY YARN 

(TWIST ADDED) 




TYPE HM 

GRAPHITE 

YARN 


300 


•‘S-GLASS 
ROVING j 
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"E -GLASS 
ROVING 


f « /C' 

' |0 iff A 

/r A 

1*1 J X* 
lui /** / *t 

U$ /,7 
*, l'* i t 



NOMEX @ 
ARAMID 


TENSILE STRAIN (%) 

ASTM 02343-67 RESIN IMPREGNATED 


STRAND TEST 

(REV 12/77) 
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the Imgth of the fibtr, «nd thlo is rsflsecsd in s loifsr «3tpsetsd strsngth 
for loQgor tsthoro. Soms dots for s limitsd rangs of L/D flbar valust ars 
shown In Flg^. 5. For Kavlar 49 » the data can be raptasantsd as 


0 


340 ( 


0 . 06 v 


0.051 


and if V tantatively extrapolate to the very long lengths contemplated. 


we calculate 

the results 

shown 

in Table 

3: 



L(m) 

0.06 

1 

10 

100 

1 Km 

10 Km 

100 Km 1000 Km 

On£(Rg/wn*) 

340 

295 

262 

233 

207 

184 

164 146 


Table 3. Detrapolated Fiber Strength for Kevlar«^9 


Clearly, the extrapolation used is queationable, and much more 
experience with long tethers is required before a firm design strength 
value can be identified. For most of the calculations in this report 
we have adopted 140 Kg/<an^ as the ultimate (break) strength, and used 
a factor of safety of 4. 

Other relevant properties of Kevlar-49 are listed in Tables 4a and 
4b. Note in particular the relatively small elongation (2.5X to break, 
or about 0.6% at the design strength used here). 

Finally, one area of some concern is the observed UV degradation 
of Kevlar-49 fibers. Here, again, the data are inadequate. Table 5 
shows a few exaixiples. The data for the 1/2*' rope Indicate partial self- 
screening, with the outer layers protecting the inner ones from the UV 
radiation. This also hints at the possibility of protective layers, 
which could also serve as a matrix for enhancing inter-fiber friction. 
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By contract, «l«ctron radiation datnag# la minimal. 
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FIGURE 5. TENSILE STRENGTH UNIFORMITY 
OF REINFORCING FIBERS 
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YARN AND ROVING OF KEVLAR* 49 
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TENSILE STRENGTH. Sj 
RESIN IMPREGNATED STRANDS 
. (ASTM D2343) 

DRY YARN 

(TWISTED TEXTILE TEST) 
VARIABILITY 

TENSILE MODULUS. Ej 

VARI^ILITY 

UNIFORMITY 


ELONGATION-TO-BREAK 


DENSITY 


VALUE 


525.000 LB/IN2 
3 620 MPa 

400.000 LB/IN^ 

2 760 MPa 

C.V. - 4 

18 X 10^ LB/IN^ 
124 000 MPa 
C.V. « 5 


2.5% 


0.052 LB/IN- 
1.44 G/CM^ 





IM 

Ii-2 

1 1-3 



1 1-4 


FILAMENT DIAMETER 
CROSS-SECTION 


SPECIFIC TENSILE 
STRENGTH. S/DENSITY 

SPECIFIC TENSILE 
MODULUS. E /DENSITY 

T 


0.00047 IN 
0.00119 CM 
ROUND 


10 X 10^ IN 
25.4 X # CM 

3.5 X 10® IN 
8.8 X 108 CM 


1 1-5 


1 1-6 


1 1-6 


TABLE 4a. 


(REV. 10/77) 


MECHANICAL PROPERTIt 









FABRIC OF KEVLAR* 49 YARN AND ROVING OF iCVLAR* 49 
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KNOT STRENGTH 



3K TENSILE STRENGTH 


FLEXURAL FATIGUE 
RESISTANCE 


CREEP. 90% ULTIMATE 
TENSILE STRENGTH 


200 CYCLES AT 56.000 
PSI OVER 3 MIL DIA. 

PIN 

(3S6 MPa; 0.08 mm) 

0.0011 IN/IN. INITIAL I 1-7 
0. SECONDARY 


COEFFICIENT OF 
FRICTION . 

YARN-YARN 

YARN-METAL 


0.46 

0.41 


FABRIC 

STRIP TENSILE 
TONGUE TEAR 
TRAPEZOIDAL TEAR 


DEPENDENT ON FABRIC 
STYLE 


II-8 


TABLE 4b. 


(REV. 10/77) 


MECHANICAL PROPERTIES 
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ULTRAVIOLET STABILITY OF "KEVLAR" A9 


1 1-9 


MATERIAL 


770 DEWIER 




(S50 FILAMENTS) 


1/8" DIAMETER 
CABLE" 


1/2" DIAMETER 
3-STRAND ROPE" 


EXPOSURE 


CONTROL 
FADEOMETER 
100 HRS 
200 HRS 

WEATHEROMETER 
100 HRS DRY 
100 HRS WET 
200 HRS DRY 
200 HRS WET 


CONTROL 
WEATHEROMETER 
100 HRS DRY 





6 MONTHS 
12 MONTHS 



22.2 

20.1 


21.3 
22.2 
18.1 

18.3 



STRENGTH LOSS 
(%) 



TABLi; 5. 

WEATHEROMETER EXPOSURE - SUNSHINE CARBON ARC 
FADEOMETER EXPOSURE - XENON LAMP 
FLORIDA EXPOSURE - HIALEAH 


^data indicates self- screening influence of 
outer layers of '*Kevlar" 49 


(REV. 11/74) 
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Appendix 2 

Platform Orbit lUcovry Uilne Iroonlelv Thruet 

Consider a platform of maae M carrying a aatelllte of mesa m, both 
of them In an orbit at (orbital apaad ). Zf the aatelllte 

la deployad on a light tether of length L» the platform descends to 


». -v^-v ' 


\ - ^ m 

p m+M 


and it travels there at 


'’l« - 


If the payload la now released, R and V become the apogee radius 

o o 

and speed for the platform in its new elliptic orbit. Since the apogee 
velocity is 


V - 

a 


2R„ 

e p 

R R +r" 
a a p 


we can now solve for the perigee radius R 


To first order in X , 

P 


we find Rp ~ 1 “ 7 Xp . The velocity at this perigee is 

■ ''« F ' \ <1 + 5 • 

In order to return the platform to its original orbit using impulsive 
thrust, we apply first a perigee impulse 


where V ' « . ^ 

" ' *p 




-Vp 


is the speed at the perigee of the new (transfer) 


orbit that will reach apogee at the intended radius R - R^^ 
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21 


Approximately, than, V' » V^(l ^ X ) , and 


* 7 V, X^ 
p 4 L p 

Whan tha platfora raaehaa apogaa at R, , a clrcularlaatlon Infmlaa 
AV. • V, - V , la naadad, whara V , - V • ^ . Wa find 

0 L A • P IV. 


AV, * 7 V, X^ 
a 4 L p 


The total AV required la cherefore 


AV - 2V. X_ - 2V, T ® 


L “’L Rj^ nrtW 


and tha total Impulse la 


MAV 


2v — 

^ L Rj^ mm 


It Is of interest to compare this impulse to that which would be 
raqulrad to place the satellite in its post-release orbit with no tether 
assist. Such an orbit has as its perigee the release radius 




M 

M*tm 


and as its perigee velocity, V - V, (1 + X ). 

P L 6 

In order to raise the satellite Impulsively from Rj^, to this 
elliptic orbit, the optimum impulsive maneuver consists of two firings, 
the first one applied at the point in the circular orbit opposite the 
eventual perigee, and such as to produce a transfer orbit tangent to the 
final orbit at that perigee (which is itself the transfer orbit apogee).* 
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Thia firing is found to bs 


% 


i \\ 


Ths sscond firing is applisd st the point of tsngsney of transfer 
and final orbits; it is found to be 


AV. 


f V, X. 
4 L a 


for a total AV of 2V. X. , and a total impulae 

L m 


mAV 


L 

IL M4« 


This is exactly the sane value found for the Inpulaa spent in 
re-establishing the platform orbit after satellite release » a 
result perhaps not unexpected on the basis of along-the-orbit overall 
oomantum conservation. 


App«adlx 3 


Platfota Orbit B» c pvry U»inn Lor Thruat 

In this nppnndlx vs «x«Bias ths orbital dynaolca of platform orbital 
rastoratlon by maans of high spaclflc Inpulsa^ low thruat anginas. Tha 
rssults will bs of usa in calculations of powar and propallant raquirs- 
manta for thasa platforms. 

1* Orbital Parturbations of tha Platferma . For tha LEO platform ralaaslng 
a payload/anglaa coubinstion, tha aaquanca of oparations ba ae follows: 

(a) Tha Shuttla docks with an orbiting platfoni which has a radial out» 
mrd tathar daployad. Payload and OTV ara tranafarrad to tha plat- 
form (including Hg for tha uppar platform). 

(b) Tha payload/anglne combination travals along tha tathar to Its top. 
Travailing rata noist ba controlled to enaure radial position at the 
and and to alnlmlsa oaclllatlons. The platform loses altitude, but, 

. to first order, the system e.g. remains in the original orbit. 

(e) The payload/anglne combination is ralassad. The platform (plus 

tathar), enter a perturbed alliotie orbit with apogee at tha release 
point. The platform mass must be sufficient to prevent reentering 
at tha perigee. 
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(d) l«ow throat anglaaa oti tha platfora ara aetivatad to alm^ly saiaa 
aad cireulatlea tha platfop orbit to ita original configuration. 


Lot ba the apogaa of tha platfora parturbad orbit* h tha tathar length 
and M and a the platform and payload/angina aaaaaa. Tha radiua of tha original 
(and eventual) platfora orbit la then 




! 


where 


X . iL. V - £ 

* R^ • M 


( 1 > 


( 2 ) 


Sene simple dynamical calculations show that, for email X * the perigee 
of the perturbed orbit ie given by 


Thus* Che eccentricity is 


,1 3 ^ 

R.+R, 2 ’ 2 l*v 

a p 


, - r- 3Xv 

R +R^ 1+V 

a p 


and the s^slmajor axis is 


R-*Wl 

*o 2 1+v^ 


or* cofflbinlug (1) * 


<3) 


(4) 


1 ^ . Xv 

-p B 1 + 4 ^ 


< 5 ) 
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A vary s^Bllar davalopaant ean ba ituid* for tko oppor (GEO) tochor nanau- 
v«r«. Tha aaquanea ia now 

(a) Tha araeion wleh a radially inward daployad tathar is initially in 
GEO orbit. The payload (af tar aaparation from the OTV last stage) 
docks with the tetbnr lower and. 

(b) Tha c.g. of: the /.yatffin is now in an elliptic orbit with apogee aome- 

what below GEO (by H, where H ie the tether length). The payload 

is made to climb along the tether at a controlled rata. At the end 
of the climb, both, the payload and the platform are in the same 
orbit occupied by the c.g. after docking (to 1st order). 

(e) low thrust anginec on the platform are activated to raise and 
eireulariae the platform-payload combination to GEO. 


It can be shown that aquations (4) and (5) still describe the perturbed 
orbit in this case, iiTith the obvious redefinitions 


X-S . v-SJESZiSSiJ 

I*CE0 “ platform) 

2. Low Thrust Steering Law . Since the action of the platform engines is 
quite gradual, we will describe their effect using the orbital perturbation 
equations (Ref . 1) . If f ^ is the applied tangential acceleration and no 
normal acceleration is applied, the rates of change of semifflajor axis a, 
aceentricity e and periapsis azimuth w are 


da 

dt 


is. 

dt 


(e cos e) 


(7) 

( 8 ) 


dw 

dt 


2f. 


ev 


sin 6 


( 9 ) 


Ref. 1 Modern Spacecraft Dynamics and Control. by M.H. Kaplan, 
i. Wiley & Sons (19>i>) , Ch. 8. 
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vh«r« 6 it agisuich from porlapalo. Tho vobiclo is a ilevly svolving 
slllpcicsl path dsseribsd instantansously by 


I _ 14^ CO S0 
r o(l-s*) 


( 10 ) 


and such that ^ 

r* II "^Maa- 0 *) <11) 

For small aceantrlclty a, we can usa (11) to allminatc time from (7), 

<8) Mid (9): 

|| « ^ ft <l-a cose) (12) 

^ ft [cos® + a(l-3coa^e)] (13) 

|| « ^ ft aiiie (1-a cose) (14) 

Tha slmplaot steering law allowing sinultanaous control of aecentrleity 
and orbital energy Is a modulated acceleration law of the form 

ft - f^ 4. cos(e-e^) (15) 

SMbstituting into (12) to (14) and averaging over one period of 6, we 
obtain (for loii^ times, neglecting products of e and or f^^) : 


da 

Om r 

de 

2a *f 

o 

y 

(16) 

da . 

a^f.cos e^ 

X 0 

(17) 

de 


dw 

SSS. a 

de 

-a 

~ f «in6 

(18) 
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For faoeoBt roduccion of oceonerlelcy, and In ordar to atraid pariapaia 
rotationt %ta cheaa " 0. l^a. (16) ,and (17) intagrata ioMdlataly to 


w U 


(19) 


•■•••".“(ris'.) 


(20) 


If at a eartaln aaloaith 6, we in^eae both a • a^ «id a « 0* we obtain the 
condition 


*1 _ "**o 


InC?^) 


( 21 ) 


So that, using equations (4) and (5) for md ^ find to 

1st o^er 


3 

2 


and aot the aecalaratibn law is 


*t “ *« a - f “»9) 

This indicates ratrofiring at perigee (f^ - i* f ) 

V A O 

forward thrust at apogee (f ^ " 4 f«) • 

i» » o 


( 22 ) 


and wucimum 
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3. Fropallant Consumption * The usual law 


*fipal 
^initial 


• a 


c 


appllaa* with 


/ ^flnal 

l*tMt 


(23) 


whara tha abaoluta value of tha applied aceelaratlon ia uaad» ainee propellent 
eonaufly>tion -la independent of thrust orientation. Eq. (11) is used again to 
alioilnate dt in favor of dS. To first order in e, the velocity increment 
par turn ia than found to be 

-1 .w. .. _..-l 2 n 



AVi - 1 2ir -4 cos * I*) + eCjV^ + 3ir -6 cos * |)3 



a 


- (7.3S1 + 7.360 e) 


(24) 


Also, to the lowest order in e» the nuniber of turns in ttee t is 

R ■ » t » so that» to that order 
Z7T * 


tv B 1.176 f t 

o 


C25) 


The product f^t can be related to the aiisslon characteristics by 
integration of the time equation (Eqs. (7.1) and (10) » combined with Eq* (19) 
for a) . Ignoring the cyclic part and retaining only the secular term, we 
obtain: 

. 3/2 


dt 

dO Hfp 4f>^2 3/4 


V^l- 6) 




^ t 1 
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which integrates to 

f,t -•4p ii - (1 - - Vf t 1 -V^ 1 (26) 


■©. O 

where (19) has been used once more. Thus, If t refers to the final elate. 
When a S| (end e 0) , we obtain 


(1 -V^ ) - *e - V 

O cL O X 


(27) 


where v is the final orbital velocity, while v would be the velocity 
®1 • ®o 

in a circular orbit with the same energy as the ixiitial (elliptic) orbit. 

a 

Using now £q. (S) for ~ , we obtain finaUy(ti7r lowest order in e) 


AV - 1.176 V 


2Xv 


(28) 


It can be seen by comparison to the results of Appendix 2 that the 
AV required with low thrust Is 1.176 times that required with the optimal 
combination of Impulsive firings. However, since the specific Impulse 
can be quite high using Ion or other electric thrustors, the propellant 
use can still be significant. 
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APPEWDIX 4 


P.vn.MMi c B of th e IcL lici ' .syrU cm du r lnr. fo wy trnrr.frr . 


^ • Po i'iPui rvob lc:u 

Wt> coriEii’oii- in liiis oct’lricn the dynamicnl effects that occur duvinr. 
arccut of a .loaded ferry vliich tranElRlct; aloiif, a tether line deployed 
frow an orbit.irf, platform. A terminal handling facility is also assumed 
to exist at the upper end of the tether; this upper platform also serves 
to provide tension for the tether, due to the ctavity gradient force 
acting on it. 

The system to be studied Is shoi<ni in Fig. 1. We vrill assume small 
angular deflections from the vertical, and ignore the mass of the tether 

itself. The latter assumption implies 
tether lengths below some 200 Km, while 
the small deflection assumption will 
be well satisfied for sufficiently small 
ferry velocity (v), provided no dynamical 
instability is encountered. These are 
precisely the issues to be clarified by 
the analysis. 

The gravity gradient forces on the three 
masses depend on orbital angular speed, 
and distance to the overall center of 
mass CM. 

With the origin of coordinates fixed at the lower platform, as shown, 
these distances are 



where 


*cm ■ + '’l*T ' y-ci. ■ '’t ^ 


X - » (l-v)x-v^x^ ; 

y M Y ^ 

^ “'em 

(l-v)y-V^L 

(2) 

*T - *«.■ 

••-’'cS 

(1-Vy) L - vy 

(3) 

^ ■ M+M 4M_ » ''t ■ 

M+M -I M 


<4a, b) 
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and 


V E 
P 


The c’aviry crnJJcnt 


] -V - 

“ 

M 

(5) 


forces 

arc 

then (po£i It Ivc upvTords) 


„r.r. 

I » 

P 


M y 
p ■'em 

(6) 

_GG 
F « 

tl! 

3R^ 

H(y -y^^) 

(7) 

T 

3R* 

Mj(l 

(8) 


In addition to these forces, the Coriolis forces must be considered, 

since the axes rotate at speed R. Wlicn the ferry is travelling upvyards 

» 

at speed y relative to the platform, since the center of mass must re- 
main (to first order) at a fixed altitude, the other masses (and the 
tether) must travel dovmwards to compensate. The absolute velocities are 
then (1 -v)y (ferry) and -vy (upper and lower platforms) . The correspond- 
ing Coriolis forces are then (positive bacl;v;ards) 

F® ■= -2R M V y (9) 

P P 

» 2R M (1 -v) y (io) 

m 


fJ - -2R Ifj, V i (11) 

We will assume the vertical accelerations are small enough that the 
tensions T^ppg^. and of the upper and lower tether segments respec- 

tively are equal to their quasi~s tat, '^c values: 


T - F®® 
upper T 


T «»-F®® 

lower p 


With the small-angle assumptions 


sin a a a: x/y ; cos a a: 1 
x» - X 

sin B = $ e — = 


L-y 


; cos B ^ 1 
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the horizontal c.oiupc>nc»its of thc*ie teiiRiiony contributo forccB 


T 

k‘ 

P 

T 

F* 

01 




V 


,CG X 
p y 


I?®® 5. 4- 

y ^ T L 
x.„- X 


T - 


L - y 


( 12 ) 

(13) 

(14) 


i'hc equations of motion for the three inassee are then 


Mp(vx + -3f)2jM[(]-v)y - V^L] + M.^.[ (1-V^)L-Vy3|^ (15) 


M( (l-V)x-VyX^.l- -3R*jM[l-v)y-V^L] + (l-VY)L-Vy]}^ 43«*M^l (l-v^)L-vyl~ 

(16) 


X™-X 


M^[(l-Vy)ji^-vxl «“3fi2M^[(l-V^.)L-vy] -ZJiM^Vy 


(17) 


It must he noticed that only too o£ these equations are independent 
since the linear combination representing the motion of the CM must be 
satisfied. The motions of ferry and upper mass relative to the lower 
platform can be extracted by the combinations (16) /H - (15) and 
(17)/M^ - (15) /Mp respectively. After simplification, 

X - ^ (vy+V^h) “ +3R® [(l-v^)L-vyl ^ - 2fly (18) 


itj. - -3R* (vy+v^L) ^ - 3P/ t(l-v^)L - yyl ^ 


"X 

y* 


(19) 


A useful variation is obtained by difference of these equations. 
Defining 


- X 


( 20 ) 


thie equation is 


1 -v. 


6 » 3R= (vy+v^I.) ~ ~ 3R* - 7 ^ 1(1-Vy)L - vyl~ 


( 21 ) 
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There nrc at t two cli.'jr;i( tcrlst:ic titie.s Involved in this problem; 
the firfJt is tlie tr.' ucit time T of. the ferry; the. Ktc.i>nd is l/fi • the in- 
verse of the orbltoJ anpulnr velocity (of the snme order a® the porJod of 
the Brav.ICy C7-ndlent OKC.Lll«tJonu) . Typically, in the- situation bcinf, 
coiiaideved, T - 1 - 3 doyti, while 1/fl 14 min. Honce* the non-dimensional 
parnmetcr 


is v<. y r.mall, and can be used as an expansion parameter for an approximate 
eolution. We make this paranete.r explicit by introducing a dimencionless 
time 



(23) 


and rewriting Eqs. (18) and (21) as 


l-v„ 




If ' V y 3 ^ I(l-vV>L-vyl —r M 


2c 




'' 0 - 3 ^ + V> 7 - 3 ^ l(l-Vj)I.-vy] i~ - Jc ^ 


V 


1-v 

V 


<36 


(24) 

(25) 


In addition to the two widely different time scales, which indicates 
the likelihood of a slowly modulated gravity gradient oscillation, Eqs. 
(24) and (25) contain the factors ~ and r“ . These will cause 

y L.-y 

singularities near the initial and final times. Physically, such 
Blngularitles arise because of the high frequency of relative oscillation'’, 
when two qf the masses come close to each other; for the condition y **- L, 
there is the possibility of a divergence of 6 as the ferry approaches the 
upper platform. 

Although Eqs. (24) and (25) are linear In (x,6), their complex struc- 

e 

ture (particularly elnce y(t) Is arbitrary) , indicates the necessity of 
approximate methods of solution. The plan of attack will be to use a WK£ 
eolution away from y»0 and y*L, and to match It asymptotically tc "inner" 
solutions valid near each end. 
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^ • The \?Lvl 1 s olut i o n (ferry not ntvir the endn) . 

Kxcepl for the very /ast osciJlations near thn endo of the trip, v»c 
onticlpato thu solution to cnnslBi of slowly modulated (on a acalc 6 • 1) 
gravity gnvltcnL nncillntlons (of period - c). There should be actually 
tv .'0 j»i>'flvlty pradii'ufc ''modec,'^ rouchly correcpondJnji to a collective., tier.r 
straight-line opcillntion of the three mnsses, and a "banding” oacl.llation 
with X opposinn and x^. The WKh method is veil suited to this liiionr 
prohlcin; v:e represent the hotTiOi»,enQOU !< approximate solutions as exponentials 
ot truncated series in c, the leading term (of order 1/c) being imagiitary 
to represent the oscillatory behavior: 


« . A(0) + EC(0) + ... ] (26) 

6 . ^ + ... 1 (27) 


where the functions K, A, B, C, D arc presumed to be smooth on the scale 
of 6. Differentiating end substituting Into Eqs. (24)^ (25) we obtain 

c* t- ^ + 21 7 (A+EC+. ..) + i|+X+etI+... + (A + CC + ...)*] - 

& & w 


- -3 y * ^ V 


(l-v^)L-vy (B-Oe... 

T=y ® ® 


(28) 


A2 K * * IT**** ** 

E*[“ fa + 21 7 (B+EDf...) +1 7 + B + eD + ... + (B + eD ...) ) - 


V 


L-y 


(29) 


where the inhomogeneous terms ± 2 e ^ have been omitted. In the under- 
standing that a particular solution will have to be added later in order 
to obtain the g«ieral solution. Here a dot Is meant to represent d/dO. 

We first observe that these two equations can be compatible to order 
only if the two right hand sides are identical (with the exp((C*’B)e) 
terms omitted). This condition leads to a second order algebraic equation 
for X se^"®: 
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/ (A-IO+ , _jr., i:\ , a-« _ _y_ ;^i . ^ (jg, 

^ V L-y vy l- v^.l. \» ^ L-y vy + V^L V * 


V 


v.lilch l»ns the Ivo ^oliirlonH 


v+ 


r. e 


X = e 


+ "^ V 

A ~B ^ 

L-y 

A'-b"" (l~Vy)L~vy 


(31) 

(32) 


If thlB coi^dltion is sntiafi^^d, the zero'th order part of, for 
instance, Eq. (29), reducett to 

V y V L-y. 

Substitution of either Eq. (31) or Eq. (32) here, leads to 
the two poseible inetantaneous frequencies: 


(33) 


- ±^f^ 


(34) 






vy+v^L 


'.j.)L-Vy] y(L-y) 


(35) 


Tlicse esqtreealons are valid for arbitrary elimbout laws. The phases 
K^, K" are obtained by time integration, and depend, therefore, on the 
particular choice of climbing law: 


k’*’ - ± 


->/r 0 

_ r r vy+v L 

-j 'V^ 7 <l^ 


de 


(36) 

(37) 


The first of these inodes is recognizable as the ordinary gravity 
gradient oscillation, at frequency'll* (2 ; we expect both X and X-, to 
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bo In pbnao Jii ili.lf! mode. Tlio cc'cond one linr: ti marc complex ntructurc, 

I 1 

with freqiicmcy Jncicaa.tnij na near y » 0 end no nenr y ■ L; 

V y ^’y 

ns v;lll be seen In the next section, there "Inner limits" of the "oittcr 
solution" lu(le«.>d uutch the outer limits of the Inner solutlono n«ar each 
exUTutc, This second ii»oile, therefore, can be expected to be the bendiiic 
tKide, and X and f should he in counterpUose. 

To continue the solution, v;c vjrltc down the order ~c parts of Eqs. 
(?6) and (29): 


2iiUi + lU « 
2iK£ + ik - 



(1-V^,)L-Vy 

L^iy 




3 


\ Vy+V^L 

V “y 




(38) 

(39) 


By division, we can eliminate (B-C) and obtain the required 
connection betvieen A and B: 


2KMK , _ ^ t-Z 2(A-B) 

2kk+K '^T <l-v^)I-vy 


(40) 


**_a A n 

For the collective (+) mode, k « 0 (aee Eq. (34)), and e^” 
ie given by Eq. (31), from which 




y(L-y) 


(41) 


dB 

aubetituting in (40), we can aolve for (aif(er cancelling dt) 

as 


dB 


(vy+Vy)L 


3? * ■'’p : (L-y r s(:y) ' 


which integrates to 


B" 


(L-y) 


^/v^(l-v, 
^ S(y) 




(43) 




<j.(i-v 

S(y) 


( 44 ) 


and, after Eq. (31), 
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This shows thst In tha collaetlva gravity gradlant ooda, tha aoplituda 
of tha oaelllatlona of X Incraaaaa naarly linaarly with y (a amplltuda about 
conatant), whlla thoaa of X^-X ^acraaaa In waplltuda about linaarly In L-y 
(6 about eonatant). Tha X-uq>lituda approachaa a finlta limit a»y 
and alnca X(y » 0) ■ 0, 6 approachaa tha aama limit aa y 0, X and X^>X 
oaeillata in phaaa. 

For tha "banding" (-) mode, for which k” is given by Eq. (35) and 

•• 

K does not v«ileh, a oomawhat more elaborate procedure la required i 
Notice that the left hand side of Eq. (40) can be written aa 


2mx 

as #• 

2KA4K 


2K^B*SK ^ 2d6 ^ ^ (2g‘Htn|K|) 

2K^A4>^ 2A+ (2A*»-An|K|) 


2A”+ 4n|K“| ■ F ; 2B"+ in|K"| - G 


(45) 


so that we have, from Eqs. (40) and (32) 


My' 'dy' v y v„L + vy 


(45) 


and, from (32) and (45) 


V (1-v )L-vy 

F - G + 2(A’ -B“) - 2 an ^ ^ 


^ ^ , -2vL 

dy ’ dy [(l-v^)L-vyl (v^L«K>y) 

dG 

Equations (4b) and (47) can be eolvad for ^ and , which can 
than be integrated to obtain F and G. a” and B*“ then follow from 
Eq. (45), using (35) for |k*~|. The result of this calculation is 
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^1/4 


(Af.) 





^ (-V) 


3 /A ■'’/Ayl/A ^L-y) 


1/4 


K3-Vi--vyj’^\/s(>') 


(49) 


wliere S(y) io as do.J!.incd In Eq. (4?). Thc.ee formulnc show that X and 

X.--X indeed osciljate In counterphase ("bending” mode). Both variables 

^ 1/4 

have aiitplltuden varying roughly as [y(L-y)) , which indicate angular 

amplitude for a like 1/y^^^ (near y 0) and for B like l/(L-y)^^^ 

(near y • L). Although this looks like a divergent behavior » matching 

to the near-end solutions will show that at least one finite angle 

Bolution exists at each end. 

A u i K/e 

Having determined e , c and e for each of the two gravity 

gradient modes, ve have a good approximation to the homogeneous solution 
and can truncate the expansion in powers of £, The r^naining task is to 
generate a particular solution of the complete equations; none can be 
identified by Inspection, and so, the siethod of variation of parameters 
must be resorted to. To this end, let us represent the homogeneous 
solution in the form 


X “ Cl fj + Cs fa + Cs fs + Cl* f^ 
^ gi + Ca 8a + C 3 g3 + Cl* gi, 


where 



- ■ 

•« e c ; 

fa 

+ K 

A -i~ 

® c e ; 

fa - '*'^1 ; 

If A 

fl* “ e 

8i 

B++1IL. 

» c € ; 

ea 

» e e ; 

g, > ; 

• 

B 

g4 “ e 


(50) 

(51) 


(52) 


( 53 ) 
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II r, , (\ , f, , f’l, ,uo » V MOW (unolinmi nf tlvi'', the 

u'l.il |msvj‘.sef iliffeu ill Jatu'ii, cnf^'rrer.ouL of aiixU l.ity c«>ndi I ions 
nn«' la'I' 1 1 1 1 a t 1 1 ' ' iu(o I ha v»>iplt'it' cciOiil ionn (I'tia. (24), (?a) .Icadti 

• 4 * * 

la .*> iVsiiM a! ai’Uvil Iou'’> I m I'l , • (•? » : 


1> (Ij I 
P,l C.\ 4 
i\ Cl 4 


f- C;« d <1 Cl 

• t 

hr C;. -I f’,.1 Cj 

• t • • 

f ? C ; 4‘ f 3 C 3 


4 fl. Cl • 0 

♦ f,ii Cl, «• 0 

4 fu Cu • ?*'- 

c< 


(5A) 


Pi C| 4 S2 Ci 4 |j3 Cs 4 iu C4 **“2^ 


B -A 

Kotico, from Kt|s. (52), (53) Hint' Bi/fi “ ga/f? “ e , and 
slwUnrly ps/fs » gti/f*. ’ • Ums, the Cirat two equations of 
(54) cdrilt tho Binplc nolutlon 


• • 

fj C|4f; Ca* 0 

• » 

f3Cs4f(,Ci,"0 


(55) 


Elimination of Ca and C«* between these equations and the bottom two 
equations of (54) then gives 


(ft - b + <a - ft’ =■ ’ 4 

'‘*1 

Noting that 

- ft’ ■ 'ft ■ ft’ •“ - ft’ **•'* 


and, eimilarly, 
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Kq5i. (.'>C') can l>c nc»l vo.J to 


A "2 . 

Cl *■ -T y 


,1 - 1/X’ 


t ' r i.n f 1 / f ? ^ ,1 i X 

T5T-)(-h. - 


f 2 • 

e ^ , ,it a,n f j/f ' V /] 1 N 

f = dO - > 


whore X , X“ arc ns In Eqs. (31), (32). Also, then 


f-h r, • fii. «• - C, 

fa 


Ca “ - Cl j C»i " - Cl, 


From Eqs. (52), (53) it follovis that ^ 


aad ^ ■" * Also, fa Is the complex conjugate of £ 

* fa “ fi* 

and, similarly, fi, “ fs . Therefore, Eqs. (57) and (59) imply 


fg Cz »' Cl 

and, similarly, Ci, » C 3 . . 


Using the expressions found before for X, K, etc., one can now 
calculate 


Cl - Ki + i.y.- f 

y3v^(l-V^)Jjj 


y ~i 1 e V ... 

e e dy 


where Ki Is an arbitrary (complex) constant and 6 is to be regarded 
for Integration as a function of y, which implies specification of 


(57) 


(58) 


(59) 


(60) 

(61) 


( 62 ) 
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A particular climbout law. Similarly, 


Cs " K 3 - iv, 




y Jk'( 0)I ^ ya-y) 


,1/4 __djr 


(63) 


V7here K”(0) is piveii by Eq. (37). 

* A 

Finally, since Cjfz “ (Cjfj) attd Ci»r«» “(Cs^s) * 


X 2 R (Cifi) + 2 R^(C3fa) 

C ” 


(64) 


2 R.(Cifj) 2 R (Cafs) 
. £- + -2- 


(65) 


or, after some reduction. 


X - 


=*’' lPco8(^e + 




4 ,) + ^ ( 6 le> 1 + 

^ Jo 


^ J COB pin 

Y S(y) 


„ .V y'(L-y») ,1/4 , _dy_l_ | 

b Q(y’)(L-Q(y')) ^ rriTT: » 


Vs(y^ 


( 66 ) 


5 / p cos (•““ e <|)) + . 

v^> 


f- 


.2 2j>i2l2l£^‘(SfJ^ll/‘lRCOB ( ) 

''t ^ ’ 


J (67) 
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Q(y) ” ^ vy 


(6S) 


and P, R, and ilf are u new net of arbitrary constuntB, to 

be found by iwpoainp, the. correct bomidary coudi Lions. 


3« Bound a ry conditions . Behavior v.’hcn the, ferry Is near one end . 

For our problem, we will nssunie the tether is initlnlly deployed 

along the redial direction, and that the ferry fitarts out from the 

lower platform vrith a relative velocity y(0) ’*= v-* in a direction 

fi-lm ^ 

making an angle (y) to the local vertical. Thus, the 

initial condltiona are 


X(0) - 0 (0) - apVQ 

8(0) ^ 0 ~ (0) —aQVp 


(68b) 


We first notice that, from Eq. (67), for 6(0) ^ 0 we need 
T cos ({> ■ 0, and since P 0 is required for later matching, we 
take 

- I (69) 


\' 3 Y 3 

such that cos (“tt” 6 + <j>) *" sin (“^ ©)• The limiting behavior 

6 C 

of the solution for small 6 is then 




2y 




Y 3 T -'0 


vT 


sin ~ (O' -9) dy’ } + 
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4. 2 


1-V, 


3/4 


,V f l\,I/4 1/4 .K (0) 

1 3 (-“■) 1 y i >1') “ 


■'T 


(i'lTjl/A 
^ 3L^ 


I 1 


3/4 




(6910 


6 « 


2L 


V *^T 

V 


-V») ^ 


h 


- 2 




E I 

4 '“3 ^ '“L^ 1 


.} (70) 


Here, the limiting form of K**(6) can be found from Eq. (37) as 




(71) 


E V3 

The integral in the first bracket is simply "-v^T (1-cos — r 0) , 
and vanishes like 0* for small G. Thus the vhole first terra in the 
expression Cor x (the collective mode) is at least of order 0^, and 
does not contribute either to x or to ^ near 6 » 0. The integral 
in the second bracket is more involved . As suggested by Eq . (71), 
we can define a time scale for the fast "bending" oscillations near 
y s 0 as 


‘'O " 3V^(l-v^) fi2L 


(72) 


Such that 


OO) 


«2vr 

•g% 


2V 


^0^0 


The integral is then 




.)lMgin (->/P-yf7 )) dy' 


^0 


(73) 
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The time T« Is much shorter than even the £ast gravity gradient 
*■ V '^0 ^ 

time • 1/n, since, from (72), 3^^ ( ' i -v T ^~L ^ C e. Therefore, 

mi intermediate, time scale exist s such that ~ » but still 0 •" |r 0, 

Wc are thus justified in cvnluatinj; (73) in nu as>nnptotJc form, for 

large values of ri » • 

''O 0 


(VoTo)^^^y ’(n’)^^^8ln(2yn’- 2</n)dy’^-(vQTp)^^V^^ 


.,^5/y/2,3/4 


(7A) 


Ve thus obtain 


* t^O ^,3/4 *3 1 ^ ^ cos (2^JT + \p) + 2 3^ (i^v-) IE ^ 


and 


n ' V Bin-Js Qt + (l-cos\3 Dt) - 

V^d-Vj.) 


Iv, 


m 1/A 1/A 2v V v't 

t Eco»(2^^+»)-3^pyrs|- (76) 


T'“ 'T 


Now, to complete the determination of constants, wr. need to 
examine the behavior on the very sliort time scale Tq ~ e^T, 
where our VJKB approximation must fail. For this purpose, we go 
back to the basic equations (24) and (25). In (24), for very short 
time after the start of the climb, we can replace 

® “ ®» 3e ® y « 0 ; 


£8 £x 3 _2L. « 2 E ^ 

® d02 ^ V v^T ^ ^ 


( 77 ) 
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The hoi;.of;('u<‘our. pnrL of (77) 1b n rcysel oqu.Ttlon In tbn v/iriabJo 
a/o . A particular solution is obtuinetl by Inspection, witJt x -> 0. 
altoijcthcr: , Uion, one obtains the solution 


X « 2fhr„T„t + D. /T 3 1 (2V “ ) + E\A' Y, (2‘/f > 

U U \ • ‘.ft 


(70) 


v,’hcre is ^ 1 S OnfineJ lu Eq. (72) „ Near the origin. 


r- 


It 


y t V. (2 V -“ ) -VTft » -F • Thus, to 

ensure x^Ostt^O, E must be Jiero, and ve have 


X • 2fhTpTQt + Ji(2',^ ) <2 QVqTjj + "p ) t 


t**0 


V^o 


(79) 


Equating the coefficient of t in tq. (79) to (See Eq. (68)) 


gives 


» - - ^0^ 


(80) 


Now, in the intermediate limit “ ® (but etill “ 0), valid 
for very Bmall c, we can use the knovm as^ptotic expanaion of the 
Bessel function Ji to obtain from Eq. (79) 


X #■ i»r Ik n ***" 

X t/tQ-» « 0 0 Y? 


(tTQ)^^^COS (2 




(81) 


This is the outer limit of the inner solution (Eq. (79)), and it 
must coincide with the inner limit (Eq. (75)) of Che outer solution. 
It can be seen that the term linear in t is already matched; matching 
of the oscillatory term requires the two conditions 




3ir 

4 


(82) 


and 


2^15? ’ 


( 83 ) 


Ill 


or, iiBiup Kqw. (f'O) I'ind (?2), 
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<jjy} (pQ- 


2 V 

3 u;jra-v^)riL, ^ 


<84) 


Only thcB constnnt P reiniilnii to be dcUrniinod , and thi© miict be 
done using the re, lining init jnl condition 

Ing iron Kq. (70) one can follov; n iuatchii\g procedure similar to 
that which led to (84). The result, after some algebra, is 


V 


V„ 0 iv 


_ W T 0"0 
2^j3 ^ 


(85) 


The final expressions for x and 6 ■> x,^~ x are therefore 

H -- y 

^2 ^ sln( jsQt) + 2 (bIhI } 

1-Vi a ' Jj, S(y} J 

« 

. v,L X 

* -F> - 


-f 


1/k 


sin t£»(K-(0-) - K-(6))I Ij ^IjyTj^ltw TT'’ 

A. H L-y / .ri 

^ \ Vi-''! i 

, rX y(w) i^'Vi \ 


- V_ 


( 86 ) 


(87) 
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^ • n i r.cn sy; lo’n r( t hi ^ ronnl L>‘. . 

ICxnminatlv'n c»f tliosc cxproMsions Cor y opproncliing L shows tlsat 

(a) the loi X cC the fcviy ripjn'oachos r.omo finite I'JmlL, with gravity 
gratlJont colU'it.fve osclllnt ionii , also of finite, amplitude. The 
''bending" rcidt/ oscillations decay as 

(b) The dif Ccrealial lag has a linoarly vanlshine component corres- 

ponding to gravity gradient oscillations v'hich tend to a constant angular 
(P) aiiplituJc. It also hfu a "bemling'' mode oscillation vrhich decays 

in ajuplitudc as (L-y)^^^; as \;e have found from the similar analysis 
neat y " 0, these constitute the asymptotic "tail" of a ncar-end behavior 
characterized by either a Jj or a Yi Bessel function. 

While near y » 0 the Yj component was absent due to the finite 
initial conditions, there is no guarantee of a similar absence near 
y •« L. A detailed analysis shows solutions of the same sort as near 
y 0, i.e. , 

5 « (i V.J.- 2nv^i,p) (T-t) F-/JT Ji + G 

where and arc the values approached by x and ^ near y * L, 

and, similar to the definition of Tq , the fast local time scale is 

' 3^ ' t wy r m 

The values of the constant F and G would now be entirely 
determined by asymptotic matching to the known solution given by 
Eqs. (86) and (87), and would therefore depend on the cllmbout law 
y(t) adopted. While the Ji part leads to oscillations of finite 
angular amplitude, the Y) part would give a finite limit for the 
amplitude of S oscillations, and hence an angular divergence. 

The object of speed control in the terminal phase of ascent should 
be to ensure the absence of this divergence. We can easily match 
(87) and (88) in their common region of validity, by noting that 
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Jl(X) ./“T COS a- ^ ) v.l,Uc Y,Ol) 

X-»o, ' * •>«> 

by cowpMlKon to K<j. (86), in onk-r t:o ensure the ab«cncc of the Yi tern, 
the phase fj'j K~(T) nu.‘«t npproauh an integer nunber of cycles: 

U T if (T) ■“ 2 IT n (50) 


^/l- flin (X - ). Thus 


ifX 


or 

V 

0 

Since this phase is varying rapidly near tlie end, fine control 
of y(t) in that phone ehould be sufficient to ensure this c edition, 
and hence to prevent uraparound. 

It nay also be noted that even the bounded fast oscillations nay 
be avoided altogethe.r near bot)t ends if the angle of departure, cIq, 
and of arrival 6^ are related in the appropriate way to the correspond- 
ing velocities and v,„. For the departure phase, this is obvious 
by inspection of Eq, (80); the condition for anooth take-off is 




^ /l du . 2 ,r n 

y V (L-y)y 


(91) 


2 a Tq 


2 ^ V ^ 

3 v^a-Vj) Qi 


(92) 


For the arrival phase, a similar sinplc expression can be arrived 
at: fast temlnal oscillations are avoided if 


^2 V , 

" t “ " 3 v^(l-v_) flL 

P p 


(93) 


5. Sone nunerical estinates . 

The order of nagnltude of the various quantities 
be appreciated by naans of a representative numerical 
the following ease: 


involved can best 
example. Consider 
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L - 200 i;m 

R *» 1.15B:;10”^ r:ul/hcr (300 Km orbit) • 

V V,^, * 0.1 , Vp 0.8 
T " 2.59;b;10^ bcc (3 Anya) 
e « ^ - 3. 332x3 O"^ 

'tloan nsccnc velocity V = 0.7716 in/sec 

— ^7 

Anc3e for smootli starting (at V) ; « 2.468x10 rad “ 

Fast Initial time scale 1' (at V) : t. ■ 1.066 sec 

Gravity gradient period - ■ •» 3133 sec "52.2 min 

/3R 


0.141' 





